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Abstract 



We investigate shock-wave solutions of the Einstein equations in the case when the 
speed of propagation is equal to the speed of light. The work extends the shock 
matching theory of Smoller and Temple to the lightlike case. After a brief introduction 
to general relativity, we introduce a previously known generalization of the second 
fundamental form by Barrabes and Israel. Then we use this to develop an extension of 
a shock matching theory, which characterizes solutions of the Einstein equations when 
the spacetime metric is only Lipschitz continuous across a hypersurface, to include the 
case when the hypersurface is lightlike. The theory also demonstrates an unexpected 
result that the matching of the generalized second fundamental form alone is not a 
sufficient condition for conservation conditions to hold across the interface. Using this 
theory we then construct a new exact solution of the Einstein equations that can be 
interpreted as an outgoing spherical shock wave that propagates at the speed of light. 
This is done by matching a Friedman Robertson Walker (FRW) metric, which is a 
geometric model for the universe, to a Tolman Oppenheimer Volkoff (TOV) metric, 
which models a static isothermal spacetime. Then our theory is used to show that 
the matched FRW, TOV metric is a solution. The pressure and density are finite on 
each side of the shock throughout the solution, the sound speeds, on each side of the 
shock, are constant and subluminous. Moreover, the pressure and density are smaller 
at the leading edge of the shock which is consistent with the Lax entropy condition in 
classical gas dynamics. However, the shock speed is greater than all the characteristic 
speeds. The solution also yields a surprising result in that the solution is not equal 
to the limit of previously known subluminous solutions as they tend to the speed of 
light. 
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Notation and Conventions 



R^ av g denotes the Riemann curvature tensor. 
Rap — R^a^p denotes the Ricci curvature tensor. 
R = g a ^R a /3 denotes the Ricci scalar. 

G, G af 3 = R a/ 3 — (1/2) g a pR denotes the Einstein curvature tensor. 

M denotes the n- or four-dimensional spacetime manifold. 

E denotes a hyper- or shock-surface in n- or four-dimensional spacetime. 

T p E denotes the tangent space of E at a point p in E. 

T P M denotes the tangent space at p in the spacetime manifold M. 

X a denotes a basis vector of T p E. 

X denotes a tangent vector in T P E. 

n denotes the vector normal to T P E 

N denotes a vector transverse to E, which cannot be in T P E. 

[ ■ ] denotes the jump in a quantity across E. For example, [g(p)] = g L (p) — g R (p) for 
p in E where g = g L on the left side of E, and g = g R on the right side of E. 
K(X) = -Vjii denotes the second fundamental form. 

K,(X) = -VxN denotes a generalized second fundamental form which depends N. 
C k denotes a function which is at least k times differentiable and it kth derivative is 
continuous. 

C k ' 1 is a C k function whose fcth derivative is Lipschitz continuous. 

diag(— 1, 1, . . . ,1) denotes an n x n diagonal matrix with the first entry equal to — 1 

and the remaining non-zero entries equal to 1. 

Notation for derivatives: For partial derivatives dF/dx a = F a . For covariant deriva- 
tives V a F = F. a . 
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NOTATION AND CONVENTIONS 



CHAPTER 1 



Introduction 

In this dissertation we give a general theory of shock matching in the lightlike 
case, and use this theory to construct a new exact solution of the Einstein equations 
that can be interpreted as an outgoing spherical shock wave that propagates at the 
speed of light. The general theory extends the shock matching theory of Smoller 



and Temple [[14], [L5| to the case of lightlike interfaces. Based on this new theory 
we construct our exact solution by matching an Friedman Robertson Walker (FRW) 
metric to a Tolman Oppenheimer Volkoff (TOV) metric across an outgoing radial, 
lightlike shock wave. In this exact solution matter crosses the interface, but nothing 
propagates at the speed of light except the shock wave. As far as we know this is the 
first such exact solution in general relativity. 

In this exact solution the shock wave emerges from the FRW origin at the instant 
of the Big Bang, and propagates all the way out to infinity. The pressure and density 
are finite on each side of the shock throughout the solution, the sound speeds, on each 
side of the shock, are constant and subluminous. Moreover, the pressure and density 
are smaller at the leading edge of the shock which is consistent with the entropy 
condition in classical gas dynamics [[/]]. However, the shock speed is greater than 
all the characteristic speeds, see [ 15|| . Subluminous shocks with this characteristic 



condition were constructed in 15 



In the general theory we translate generalized notion of the second fundamen- 
tal form of Barrabes and Israel JTJ into the shock matching framework of Smoller 
and Temple [13, TE\. We base the analysis on a modified Gaussian Skew (MGS) 



coordinate system on lightlike surfaces which we introduce in this dissertation. 
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2 1. INTRODUCTION 

One surprise is that the relation between the sound speeds on the front and back 
sides of the shock in our new exact solution, are not equal to the speeds in the limit as 
the shock speed in the Smoller- Temple solutions, obtained in [ 15[| , tend to the speed 
of light. 

Another surprise is that the general theory shows there exist gravitational metric 
components which are C 1,1 across a lightlike surface for which divG = does not 
hold in the weak sense. This implies that the matching of the (generalized) second 
fundamental form alone is not a sufficient condition for conservation conditions to 
hold across the interface. 

1.1. The Central Problem: Lightlike Shock- Waves 

The goal of the work here is to construct shock-wave solutions of the Einstein 
equations which move at the speed of light. In Einstein's general relativity all physi- 
cally possible spacetimes correspond to solutions of Einstein's Equations 

(1.1) G = kT } 

which are metrics that describe the geometric structure of spacetime. Equation fll.lQ 



is a tensor equation representing 10 nonlinear partial differential equations where G 
describes the geometric structure of spacetime, and T represents the matter which is 
the source of the gravitational field. 

Shock-waves were first studied in compressible, non-viscous gas flow as discontinu- 
ities that form in the fluid quantities pressure, density, etc. The mathematical theory 
of shock-waves is contained in the study of hyperbolic conservation laws, and applies 
to much more that just gas dynamics. Mathematical shock wave theory not only 
models phenomena such as the sonic boom created by a fast moving plane, but the 
same theory also incorporates the propagation of "gridlock" in traffic flow, the leading 
edge of a nuclear explosion, flame fronts in combustion, and separation of boundaries 
between chemical species in chromatography ||17|| . If we take the covariant divergence 
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of equation flLip , then it turns out that 

(1.2) divT = 0. 

In the limit of low velocities and weak gravitational fields equation ( |1.2|) reduces to 
the classical compressible Euler equations which is a hyperbolic conservation law in 
gas dynamics. 

The shock-wave solutions considered in this dissertation are constructed by match- 
ing two metrics (geometries) across a surface embedded in spacetime. In their well- 
known 1939 paper | |1 1|| describing gravitational collapse of a star, Oppenheimer and 
Snyder gave the first example, which had interesting dynamics, of a solution of the 



Einstein equations obtained by matching two solutions across a surface ||14j| . With 
their simplifying assumption that the pressure is zero, the surface across which the 
metrics are matched is not shock-wave, but a contact discontinuity which is a discon- 
tinuous solution where neither mass or momentum cross ||14| . Smoller and Temple 
extended the Oppenheimer-Snyder model of gravitational collapse to the case of non- 
zero pressure in their 1994 paper [[L4| . Their theory was based on work done by Israel 
in his 1966 paper ||], which related the second fundamental form across a shock 
surface to the Rankine-Hugoniot jump conditions across the shock. Then, in later 
work ||15|| , Smoller and Temple constructed a family of exact, spherically symmetric, 



shock-wave solutions of the Einstein equations by matching two spherically symmetric 
metric across a surface. The shock surfaces in Smoller and Temple's work are assumed 
to be non-lightlike, that is, moving slower than the speed of light. The contribution 



here will be to extend their results given in [ |14|j and ||15|| to also incorporate the case 
when the surface is lightlike. 



The Two Difficulties of the Lightlike Case. In the lightlike case the mathe- 
matical machinery used in the sub-luminal case breaks down in two areas. First, the 
induced metric on a lightlike surface is degenerate, and the second is that the second 
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fundamental form, 
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(1.3) K = -V x n, 

cannot be used to describe the dynamics of the surface in the ambient spacetime. A 
lightlike surface is characterized by the length of the normal vector being zero. Note 
the in this case the normal vector is non-zero, but the metric is not positive definite. 
A normal of zero length is orthogonal to itself, hence it lies in the tangent space 
of the surface. Since the normal is no longer transverse to the surface, the second 
fundamental form no longer gives information about how the surface in embedded in 
the spacetime manifold. 

The degeneracy problem is dealt with by considering the problem in the context 
of the whole spacetime manifold where the metric is not degenerate. Unfortunately, 
the failure of the second form is not fixed so easily. To rectify this failure we make use 
of an idea originated by Barrabes and Israel in Jl| in which they define a generalized 
second fundamental form in terms of a vector N transverse to the surface. 

The work of Barrabes and Israel in |]J focuses on the dynamics of surface layers 
whose theory is also based on matching two metrics across a surface. A surface layer 
differs from a shock- wave in that a shock surface is characterized by a jump in the 
density of the fluid across the surface, in contrast, the density becomes infinite in a 
surface layer. Mathematically, the Einstein tensor G, which comprises the right hand 
side of equation ( |1 . 1| ) , contains no delta-function singularities across a shock surface, 
but across a surface layer G does contain a delta-function singularity. We also note 
here that Barrabes and Israel use a scalar version of the second fundamental form 
where we use the form given in equation (|1.3|) which is a tangent vector in the surface. 



Regarding the lightlike case, Barrabes and Israel state that because of the breakdown 
in the second fundamental form, the lightlike case is a relatively neglected area of 
surface dynamics which remains imperfectly understood ffl. 
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1.2. How This Document is Organized 

In chapter |2| we give an introduction to general relativity, and a derivation of the 
second fundamental form. Then in chapter ^ we give a generalization of the second 
fundamental form based on the idea of using a transverse vector in place of the normal 
vector. Continuing on, we state and prove the main result of this dissertation giving a 
set of equivalent conditions which yield the existence of shock solutions in the Einstein 
equations when the metric is Lipschitz continuous across a hypersurface, and we finish 
the chapter by proving a similar result involving spherically symmetric metrics. In 
chapter § we give an exact, spherically symmetric, lightlike shock-wave solution of the 
Einstein equations based on Smoller and Temple's work in | |15[ |. Chapter [5] concludes 
the results. 



INTRODUCTION 



CHAPTER 2 



General Relativity and the Second Fundamental Form 

The purpose of this chapter is to make this dissertation self-contained for anyone 
with a good understanding of differential geometry. We construct the geometry of 
spacetime from physical assumptions, derive the stress-energy tensor for a perfect 
fluid, define the Einstein equations, and conclude the chapter with an introduction 
of the second fundamental form. 

The underlying theme of general relativity is the idea that the physics of our 
universe and its geometry are interdependent in such a way that they cannot be 
separated from each other. The Einstein field equations, 

G = kT, 

is a mathematical formulation of this idea with respect to the physics of gravity. The 
Einstein equations match G, a second order differential operator on the metric that 
is related to the Riemannian curvature tensor, to T, the stress-energy tensor of the 
fluid or matter in the region of spacetime that is being considered. 

2.1. The Geometry of Spacetime 

A theory of gravitation begins with a notion of space and time. Spacetime is 
a continuum of events each of which has three spacial components, and one time 
component. If we give each event a name in terms of these components, say x = 
(lX . Ob ■ Ob ■ Ob 3 ) where x° = ct with c the speed of light in a vacuum, t denoting time, 
and x a for a = 1, 2, 3 representing the spacial components, then we can represent each 
event in as a coordinate x. The geometry of spacetime is the structure of how this 
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continuum of events fit together. We shall see, it is the principles from physics that 
will determine this structure. 

Mathematically spacetime is a four dimensional manifold M equipped with a 
metric g. It is the metric that carries the information about the rates at which 
clocks run, and the distances between points. It is also the metric that describes the 
geometry. The question now is how do we determine a metric which models spacetime 
or subset of it? 



Special Relativity. In order to proceed we need the notion of an inertial coordi- 
nate frame or sometimes called an inertial observer. A coordinate system is inertial if 
any two events in the coordinate system share the same clock or time coordinate, and 



the geometry of space at any constant time is Euclidean ||12| . Consider a spacetime 
without gravity, that is, a spacetime containing negligible mass. This is the realm of 
special relativity which can be deduced by the two following postulates: 



Postulate 1 (Principle of Relativity) . Let x be an inertial coordinate system. 
Then any other coordinate system x which moves uniformly, and is non-rotating 
relative to x is also an inertial coordinate system. Furthermore, "The laws of nature 
are in concordance for all inertial systems" B] . 



Postulate 2 (Propagation of light). The speed of light in a vacuum c is ob- 
served to be the same from any inertial coordinate system. 



The principle of relativity is equivalent to the statement that any body in uni- 
form motion remains in that state unless acted upon by an external force ||12|| . The 
universality of the speed of light means that two observers, one moving with non-zero 
uniform velocity with respect to the other, will each observe a single light ray moving 
through vacuum to be travelling at same speed c. 
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Measuring Distance in a Special Relativistic Spacetime. Consider two inertial 
frames O and O, each having only one spatial dimension for simplicity, with co- 
ordinates (x, ct), and (x, ct) respectively. Assume O is moving along the x axis with 
relative uniform velocity v see figure [I]. Let two events be connected by a light ray. 





light ray 






light ray 


ct 




ct 












V 




O 


X 





X 



Figure 1. 

Then in the O frame the squared distance between the events is given by 

(2.1) As 2 = -(Act) 2 + (Aa;) 2 = 0. 

By the propagation of light law the same two events, in the O frame , also satisfy 

(2.2) As 2 = -(Act) 2 + (Ax) 2 = 0. 

Now, assume that coordinate transformation from O to O is linear, and that their 
origins coincide. Then if follows that 

As 2 = As 2 

for any two events in spacetime for any inertial coordinate frames O and O whose 
origins coincide, and are related by a linear transformation. 

Recall that the metric g is required to give data on the rates at which clocks run, 
and the distances between spatial points. Since any data determined by the metric 
also needs to satisfy postulate [3], we require that g be a coordinate independent 
quantity. Therefore, g should incorporate the coordinate invariant quantity As 2 , 
which composed of squared displacements between events. 
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What is the formula for g in a given inertial frame? Let A denote a displacement 
vector between two events in some inertial coordinate system. Then the squared 
distance between the two events in the coordinates is given by 



g(A, A) = -(A ) 2 + (A 1 ) 2 + {Ay + (A 6 ) 



1\2 



i2\2 



t 3\2 



which consistent with the squared distance in fl2.2|) . This leads to the inner product 
between any two vectors A, and B given by 



g(A, B) = -A U 5 U + A B + A 2 B 2 + A 6 B 6 . 

Notice this is very similar to the standard dot product with the only difference being 
the negative sign on the product of the first components. We can also write this 
metric in the following form: 

T 



(2.3) 



g(A,B) 



A 1 

A 2 



/b°\ 

B 1 
B 2 

3 



f-i o o o\ 

10 
10 
\ I J \B A J 

In this form we can see that g has a signature of +2. A metric with this signature 
is said to have a Lorentzian signature. We will only deal with metrics of Lorentzian 
signature. Furthermore, the metric given in ( |2.3| ) is said to be Minkowskian, and is 
denoted by r] a/3 = diag(— 1, 1, 1, 1). The inner product g is not positive definite. To 
see this consider 

g(A, A) = ((1, 1, 0, 0), (1, 1, 0, 0)) = -1 + 1 = 0. 

However, we can always find a coordinate basis so that g is of the form r] a p. Therefore 
g is non- degenerate, that is, g(A, B) = for all vectors B if and only if A = 0. 

Definition 2.1. A vector Ais said to be spacelike if g(A, A) > 0, is lightlike if 
g(A, A) = 0, is timelike if g(A, A) < 0. 
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From this we can give the causal structure of spacetime, that is the causal rela- 
tionship of an event to other events. The future or past of an event p lies inside the 
light cone of p, see figure |[ If an event p lies at the origin of an inertial frame, then 
the future and past of p is given by all events q such that the displacement vectors 
between p and q are timelike. 

Future 

Future light cone 



Timelike related 

Past 

Figure 2. 

Summary 2.2 (Special Relativity). From postulates [I] and ^| we were able to de- 
scribe how the squared distance As 2 should be measured, and that it is invariant 
under coordinate transformations between inertial frames. From the idea of As 2 we 
defined the inner product g given in Q2.3| ), and designated it as our metric for space- 
time with no gravity. From this we ascertained that g is a metric with Lorentzian 
signature. 

General Relativity. We now consider spacetime with gravity. We are motivated 
by the fact that objects under the influence of gravity move along free fall paths. To 
see this imagine an observer in a spaceship with the engines turned off, and orbiting 
the earth under the influence of only earth's gravity at a distance r above the earth. 
We are assume the earth is a perfect uniform sphere, and ignoring the fact that the 
earth is rotating. Then an observer inside the spacecraft now feels weightless, and 
can position themselves in such way that they remain stationary in the middle of the 




Spacelike related 
Past lightcone 
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ship so that they are not moving relative to the walls of the ship. The ship, and the 
observer are moving along a free fall path of earth's gravity. Such a coordinate frame 
is said to be comoving with the metric. 

How does the observer know that they are under the influence of a gravitational 
force? Consider a second ship in orbit the same distance r above earth, but positioned 
a distance d to the right. Also, assume the second ship is moving with the same 
velocity parallel to the original ship at some instant of time say to, see figure |3|. As 




Moving parallel at time t 

Figure 3. 



seen in figure |3|, the distances between the two ships is not constant, and their paths 
will even cross. In this space Euclid's parallel axiom does hold, and spacetime is 
curved. 

Mathematically, for the local coordinate system, (inside the original ship) the 
Minkowskian metric, g a p = r) a p, is appropriate just as it was in special relativity. 
Also, at the point of spacetime where the ship is located, we have g a p a = 0. This 
defines a locally inertial frame or also called locally Lorentzian. 

A logical question to ask at this point is, does there exist a global inertial frame 
when gravity is present? The answer is no. We have already seen from our example 
that the two ships cannot be described by the same inertial frame because they lie 
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in a curved space which is not Euclidean in it's spatial components violating the 
definition of an inertial frame. This violation is due to the non-uniformity of the 
earth's gravitational field. As we expect with a curved space, we can only consider 
local frames diffeomorphic to M 4 . 

The idea that objects moving in a gravitational field are described by comoving 
local inertial systems takes into account the Equivalence Principle. The Equivalence 
Principle says that all bodies move the same way in a gravitational field ||18||. It is 



the equivalence principle that motivated Einstein to formulate the theory of general 
relativity. 

The implication that under the influence of gravity particles move in a curved 
space, which can only be described in the context of locally inertial frames, is that 
we have to use the mathematical machinery of differential geometry to characterize 
the dynamics of particles under the influence of gravity. This is in contrast to special 
relativity where the curvature is zero, since their is no acceleration of particles due 
to gravity, and one coordinate frame is sufficient to describe spacetime. Therefore, 
we will present the concepts from differential geometry that are needed to describe a 
curved spacetime. 

The Metric in a Curved Spacetime. The metric completely defines the curva- 
ture of the spacetime manifold, and so we begin by constructing the spacetime metric 
g for a possibly curved spacetime. Earlier we used displacement vectors to construct 
the Minkowski metric in the flat spacetime of special relativity. If the displacements 
are infinitesimally small, then we can associate them with tangent vectors, which are 
defined by directional derivatives. Let p be in a point in the spacetime manifold M, 
and suppose 

x = 0°, x 1 , x 3 , x 3 ) : M — ► R 4 , 

are coordinates on M whose origin is at p, where x° = ct, c = 1, and x % for i = 1, 2, 3 
denote the spatial coordinates. Also let T denote the set of all C°° functions that 
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take smooth curves on M into R. 



(2.4) 



M 



Then a tangent vector X at p in M is defined to be a map X : T — > R. For any 
coordinate system a; on a neighborhood of p, there exist a coordinate basis given by 



(2-5) 



X a (/(P)) 



(/OX" 1 ) 



x{p) 



See diagram in (|2.4|) . That is, {d/dx a } gives a basis for T p M, the tangent space of 
M at p. 

We now introduce the Einstein summation convention which says that for any 
expression, equal up and down indices are summed over all possible values the index 
can take. For example, using (p.5[) , if X 6 T P M, then 

3 



x = J2 xa 



o=0 



_d_ 

dx c 



X c 



,_d_ 

dx c 



Notice that the upper index in the denominator is considered a down index. 

With tangent vectors defined as directional derivatives, they are in one-to-one cor- 
respondence with displacement vectors. Since the metric is quadratic in displacements 
we define g as a map where 

g : T P M x T P M — ► R. 

Also, we assume that g is symmetric, and non-degenerate as was the case in special 
relativity. 

The metric g being symmetric guarantees, that in coordinates, the components 
of the metric g a @ make up a symmetric matrix, and therefore there will always exist 
a linear transformation to another coordinate system that will take g a p to r] a p = 
diag(— 1, 1, 1, 1), the Minkowski metric. Furthermore, it can also be shown that the 
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metric g is locally Lorentzian or locally inertial, that is, given p in M, there exist 
coordinates x whose origin is at p such that 

9ap(p)=Vap for a, (3 = 0, 1,2,3, 
9ai3,»(p) = Q for a,/3, 7 = 0,1,2,3, 

and 

at the point p for at least some values of a, f3, [i and v if spacetime is not flat, that is, 
there exists a gravitational field. This corresponds exactly with the notion of locally 
inertial frames given in the previous section. 

Tensors in Spacetime. An alternative version of the Equivalence Principle is 

the General Covariance Principle which motivates the use of tensors to measure 
physical quantities that depend linearly on displacements. The General Covariance 
Principle states that an equation holds in a gravitational field if the two following 



conditions are met 20 



1. The equation holds in the absence of a gravitational field; that is, it is consistent 
with postulates [l] and |[ 

2. The equation is covariant, which means that the equation holds under any 
coordinate transformation x — > y. 

We take the following definition of a tensor from |18|| . Given a finite vector space 



V, and denoting its dual space by V*, a tensor, T, of type (k, I) on the space V is a 
multilinear map 

(2.6) T : (V*) k xV l — ► R. 

Note that our definition of tensor makes no mention vector or dual vector components. 
A tensor gives the same real number for a particular set of vectors and dual vectors 
independent of the coordinates the components are computed in. 
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Our first example of a tensor is the metric g which is a (0, 2) tensor, that can be 
written in terms of coordinate basis one-forms as 



(2.7) ds 2 = g a/3 dx a ® daP = g a/3 dx a dx /3 . 

It is customary to drop the exterior product sign between the one-forms dx a . Further- 
more, the metric provides a mapping between vectors and one-forms (dual vectors) 
at every point. Thus, given a vector field X(p), there is a unique one-form field given 
by X{p) = g(X(p), ). 

Using this map between vectors and one-forms we can construct a unique and 
useful one-form basis. Suppose {d/dx a } is a coordinate basis for T P M. Then {dx a = 
g(d/dx a , )} gives a one-form basis. 

Since g is non-degenerate, there exists an inverse which takes one-forms to vectors. 
If the indices of the basis one-forms are "up" and "down" for basis vectors, then we 
can use the Einstein summation convection to keep track of these mappings. Notice 
that for the Einstein summation convention to work we have to have the indices of 
vector components "up," and the indices on components of the one-forms "down." 



For example, see equation (|2.7| ). The components for the inverse of g are denoted 
by g af3 . Notice that g aa g al3 = 5f as we would expect from matrix multiplication. 
Using the metric g as an invertible map from vectors to one-forms can be executed 
by raising and lowering an index of a tensor by contracting the index with metric. 
For example, we can map the vector X to a one-form as 

ga<jX a = X a . 

Another example mapping a (3, 1) tensor to a (2, 2) tensor: 

o^R 01 = R a P 
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The components of a (k, I) tensor given in one coordinate system can be written 
in terms of another set of coordinate system via the tensor transformation law: 
(os.) rpai,...,a h _ rpni^. ^y UJb ux 

Here we have written a ^-coordinate tensor component T^'"''g k in terms of x-coordinates. 
Also, note that the Jacobian matrix satisfies || = (||) 1 . 

Parallel Transport and the Derivative Operator. Motivated by the idea of 
describing the curvature of spacetime intrinsically, as opposed to describing spacetime 
as an embedding in some other space, we will define curvature in terms of parallel 
transport. Intuitively, a vector field Y defined on every point along a curve is said 
to be parallel transported along the curve if the vectors of Y are parallel, and are of 
equal length at infinitesimally close points. Mathematically, in a locally Lorentzian 
frame at a point p, the components of the Y vectors must stay constant along the 
curve near p. If we let be a parameterization of the said curve, and denote the 
its tangent by X = dx/d^, then dY a /d^ = at p. However, 

(29) dY 1 _dxdY 1 _ dY 1 _ 

1 ' d£ ~ d£ dx dxP ~ 



Equation (|2.9| ) leads to the idea that in order to define parallel transport for any 

coordinate frame, not just an locally inertial system, we require a notion of how to 

take derivatives of vector fields. 

Consider the vector field Y in an arbitrary coordinate frame, not necessarily 

Lorentzian. We differentiate Y as follows: 

d ( d \ dY a d d d 

K ' J 13 dxP\dx a ) dx? dx a dx? dx a ' 

Now, 

d d 



dx 13 dx a 

is a vector and can be written in terms of the coordinate basis, that is, 

(2.11) AJ-=i^ ( 9 



dx?dx a ^dx^' 
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where T^p is called a Christoffel symbol which is yet to be determined. Therefore, we 
can write equation fl2.10|) as 

dY a d d d d 

(2.i2) v,yEK r -- + yT;-Er„- + yT; 



dx? dx a a/3 dx» fi dx a aP dx»' 

The differential operator V is called the covariant derivative. Given arbitrary vector 
fields X and Y the covariant derivative is defined as 

d 

(2.13) V X Y = X a V a Y = X a Y p . 



Notice that our definition of covariant derivatives did not involve the metric. 
However, if we recall that the metric g maps vectors into one-forms, then it would 
seem that the metric would have something to do with how their derivatives are 
related. In a locally Lorentzian coordinate frame we know that for any vector field Y, 

Y =Y 

since the derivatives of the basis vectors are zero just as in special relativity. The 
same relation holds true for any tensor in a Lorentzian frame including the metric. 
Therefore, in a locally Lorentzian frame, the covariant derivative of metric components 
are given by 

(2.14) 9aP;i ~ 9a@,j = 0- 

Invoking the General Covariance Principle, we find that 

(2.15) g a p. n = 

holds in any coordinate frame. Using this result any Christoffel symbol can be written 
in terms of the metric. To do this we must first show that Y a = T° . 

Consider a scalar field in a Lorentzian coordinate frame. Then V0 is a one-form 
with components )( g. In the Lorentzian frame 



(2.16) 
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since partial derivatives commute. Once more, invoking the General Covariance Prin- 
ciple, we have that the symmetry in equation ( |2.16|) holds in any coordinate system. 
Therefore, in any coordinate frame 

(f>,a;/3 - <P,f3;a = (f>,a,/3 ~ ^a^,* ~ + = 0, 

and we have that 

(2-17) K p = T^ a . 

Then, since g a /3- n = in any coordinate system, we can write 

9ctf3,~{ — T a ^g a p + r^^o-aj 
9 (3^, a — r / g a 5'o-7 + ^ la 9(r(3- 

Thus 

(2.19) — g a p tl + g ia ,p + 9p-y,a 



Using the symmetry property ( 2.17 ), equation ([2.1 9|) the Christoffel symbol can be 
given in terms of the metric as 

(2-20) Tip = -g ai {-g^ + g ia ,p + g Pl , a }. 

Now that we have a differential operator in hand we can define parallel transport 
without reference to a coordinate system. This is accomplished by expressing the 



parallel transport condition given in equation ( |2.9| ) in a covariant way. Rewriting ( |2.9| ) 

as 



(2.21) 



X p Y a j = X p Y a p = 0, 
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which holds, if and only if, 

(2.22) v x r = o, 



see equation ( p.!3| ). Since equation ( |2.22| ) is a tensor equation, it must hold for 



any coordinate system. Equation ( [2.22|) defines parallel transport of a vector Y 
along a curve with tangent vector X. Furthermore, given any initial value of Y, 



equation ( 2.22 ) has a unique solution, and so given any initial vector Y at a point on 
a curve this defines a unique vector at any other point on the curve. Then we can 
use the structure of parallel transport to map T P M to T q M along a given curve from 
point p to point q. The covariant derivative, which gives this notion of identifying 
tangent spaces from different points along a curve, is also referred to as a connection, 
and the Christoffel symbols F^g are referred to as connection coefficients. 

With a definition of parallel transport we are in a position to describe, mathe- 
matically, the "free fall" paths in a gravitational field. These paths correspond to the 
geodesies in the spacetime manifold. The geodesies of the flat spacetime in special 
relativity are straight lines. Qualitatively, geodesies can be thought of as the straight- 
est lines possible in a curved manifold, and they also correspond to the "straight" 
lines at the origin in a Lorentzian coordinate frame. In flat manifold the straight lines 
are the only curve in which the tangent vector parallel transports itself. Generally, a 
geodesic is given by the equation 

(2.23) V X X = 0, 

which finds the curves where the tangent vectors parallel transport themselves. If we 
let £ be a parameter which gets mapped to curve, then 

_ dx 13 d d 



and we can write the geodesic equation ( p. 23 



as 



d ( dx a \ dx a dx 13 

di\dj) + ct/3 aTaT 



( 2 -24) -A — )+r^— — = 0, 
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which is a second order quasi-linear differential equation for x a (£). Equation fl2.24|) 



has a unique solution given an initial position (£o)j and an initial direction 

Xq = X a (^ ) at £ This idea will be used in chapter [| that is, we will be given 



a point in spacetime and vector at that the point, and be able to define a unique 
geodesic through the point. 

The Curvature Tensor. Having defined parallel transport we are now able to 
define curvature. The Riemann Curvature Tensor describes the failure of initially 
parallel geodesies to remain parallel ||18|| . This corresponds to the situation involving 
the spacecraft orbiting earth in figure |3|. 

If we parallel transport a vector X around a coordinate grid loop starting and 
ending at point A in figure (01), then it can be shown, using the general coordinates 




x = a +Aa 



Figure 4. A piece of the coordinate grid. 

x a and x 13 , the change of X 7 by parallel transport along Aad/dx a , then Abd/dx 11 
then —Aad/dx a , and back to A along Ahd/dx 13 is given by 

AX~< = AaAb(r R -T\ + T\T° - FL ) X", 

V M«,p MPiCf crp IJ.ce era tip ) ' 



see |fL2 , section 6.5]. From this comes the definition of the Riemann curvature tensor: 



11 avf3 — L af3,v 1 av,\3 1 au L a/3 1 a/3 L aw 
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The Riemann curvature tensor R^ au p transforms as tensor, via the tensor transfor- 



mation law ( |2.8|) , giving a covariant measure of the second derivative of the metric 
9- 

Proposition 2.3. The Riemann curvature tensor satisfies the following identi- 
ties: 



(2.26) 


rL avf3 




(2.27) 


U [au/3] 




(2.28) 


R/iav(3 


= 9/mtR av/3 


(2.29) 


R/iav(3 


Rvfifjka. 



The Ricci tensor is a contraction of the Riemann tensor, and is given by 
(2.30) R af3 = R u av p. 



Due to the symmetries given in proposition |2]^ any contraction of the Riemann tensor 
reduces to ±R a p, and R a p = Rp a . Similarly, the Ricci scalar is given by 

(2.31) R = g a ?R aP . 

Einstein Tensor. The Einstein tensor, which comprises the left hand side of the 
Einstein equations G = kT, is the simplest (0, 2) tensor constructed from the Riemann 
tensor R^g and the metric g a p such that 

VG = divG = 0. 

We shall see in section 12. 21 that divT = 0, so that when G = kT we must have divG = 



0. To derive the Einstein tensor we consider the Riemann curvature components R^ ai/ p 
in a Lorentzian coordinate frame. In a Lorentzian frame the Riemann tensor, given 



in equation (|2.25|) , can be written as 

(2-32) R^au/3 = 2^(9(7(3,^ — 9av,a/3 + 9av,a/3 ~ 9af3,c 
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Lowering the index pL in ( |2.25| ), and taking the partial derivative with respect to x 
yields 
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A 



(2.33) 



Exploiting the symmetries of the metric g a p = g/3 a , and the fact that partial deriva- 



tives commute, equation (|2.33| ) can be written as 



(2.34) 



R^av0,\ + Rfj,a\v,f3 + R^a/3X,u — 0. 



Now, since our coordinate frame is locally Lorentzian we have T^g = 0, and so 
equation ( |2.34| ) is equivalent to 



(2.35) 



Rriavl3:\ + Raa\u\0 + Rn,a0\\v ~ 0. 



The relations of the components given in equation (|2.35|) are called the Bianchi Iden- 
tities, and since ( |2.35| ) is a tensor equation it is valid in any coordinate system. 
The covariant derivative is commutative with respect to contraction, that is, 



for any tensor T. Also, by equation ( |2.15| ), g a p- n = 0, and since the inverse components 
g a P a functions of the metric components g a p it follows that 



g aP ;i = o. 



Therefore, if we apply the Ricci contraction |2.30| to the Bianchi identities (|2.35|) , then 



(2.36) 



g^ V {R^avf3;\ + Rfj,a\v;l3 + R, 



p,af3\;u) — Ra/3;\ ~ Ra\;/3 + Ra0X;u ~ u ' 



0. 



where we have used the antisymmetric property of the Riemann curvature tensor 
given in equation fl2.26j ). Now, we contract again on the indices a and /3, yielding 



(2.37) 



g 0/3 {Raf3;\ — Ra\;f3 + Ra/3\;i>) ~ R;^ ~ R\;0 ~ R\;v ~ 0> 
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or 



(2.38) (2i£ - 8lR). iV = 0. 



Then we define the Einstein curvature tensor as 



^2.39) G aP = R aP - \g aP R = G^, 



and by equation ( |2.38| ) we have 

G aP . p = 0. 

2.2. Describing Matter in Spacetime 

In this section we will derive the stress-energy tensor. It is this tensor, T, which 
describes matter, that will be inexorably linked to the geometry of spacetime via the 
Einstein field equations G = kT . We assume that physical objects in the region of 
spacetime under consideration are regarded as a fluid, that is a continuous distribution 
of matter. In particular, we assume the distribution of matter is approximated by a 
perfect fluid which is described as a fluid with no heat conduction or viscosity. No heat 
conduction implies that energy can flow only if particles can flow, and no viscosity 
means that all forces are perpendicular to the interface between particles ■ There 
are other stress-energy tensors for fluids not incorporating the assumptions of a perfect 
fluid, but the tensor derived below will be that for a perfect fluid. 

The Stress-Energy Tensor. We begin with the notion of the four-momentum. 
A curve passing through an event is classified as timelike, lightlike, or spacelike ac- 
cording to whether the inner product of its tangent vectors g a pX a X^ is timelike, 
lightlike, or spacelike, see definition |2.1| . The path of any material particle passing 
through an event p in spacetime must lie inside the lightcone of p, see figure |2], oth- 
erwise an observer at p would see the particle moving faster than the speed of light. 
Thus any moving particle must lie on a timelike curve. The timelike path of a moving 
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particle may be parameterized by the proper time r which is defined by 
(2.40) t= j \-g a pX a X^f' 2 df, 

where £ is any arbitrary parameterization of the path, and X a is the tangent vector 



to the curve in this parameterization [|18 



The four-velocity is defined by the tangent vector u a of a timelike curve param- 
eterized by proper time r. In the presence of gravity any material particle has an 
associated rest mass m, which shows up in the equations of motion ||18|| . The four- 
momentum vector, p a of a particle with rest mass m is given by 

dx a 

(2.41) p a = m = mu a . 

ar 

The Energy of a particle with four-velocity v@ is given by 

(2.42) E = -g a pp a v^ 



taken from |18|. In a Lorentzian frame moving with the particle, called a rest frame, 
we have dr = dx°, and u° = 1, and u a = for a = 1, 2, 3, and so E = —gooP°u° = 
p° = m. We note that this is the well known equation E = mc 2 in our units with 
the speed of light c = 1. This We identify the energy of a particle with the 0th 
component of its four-momentum. Furthermore, this also shows that, in general, we 
cannot differentiate between mass and energy. 



Now, we define the general stress-energy tensor, in terms of its components, as 
(2.43) T a/3 = 



The rate at which the ath component of 
momentum crosses a surface of constant . 



From this definition, along with our assumption of no heat conduction, and no vis- 
cosity we can derive the components of the stress-energy tensor in a Lorentzian frame 
moving with the same velocity as the fluid. 
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• T 00 is the rate at which energy crosses the surface of constant time, and is 
referred to as the mass/energy density or simply the energy density, and is 
denoted by Greek letter p. 

• T a0 = T 0a = is a consequence of no heat conduction. The rate of momentum 
in the direction x a across a surface constant in time is zero because without 
heat flow the momentum of a particle cannot change without moving, which 
means T a0 = 0. Since the coordinate frame is moving with zero velocity with 
respect to the fluid, with no heat conduction the energy is constant across a 
surface of constant x a which means T 0a = 0. 

• T ab is the rate of the ath momentum component across a surface of constant 
x a . The absence of viscosity, a force parallel to the interface between particles, 
implies that the force should be perpendicular to the interface. Consequently, 
a surface of constant x a will only be the force per unit area against it in the 
x a direction which is equal for each a. Mathematically, this means T ab = p5 ab 
where p denotes the pressure. 

Therefore, in a Lorentzian frame comoving with a perfect fluid the stress-tensor is 
given by 

(2.44) T a(B = pg a/3 +{p + p)u a u p , 
and written in covariant form 

(2.45) T = (p + p)u®u + pg. 

Conservation of Energy and Momentum. The stress-energy tensor T is a repre- 
sentation of the energy and momentum in a fluid, therefore it seems reasonable to 
expect there to be a way to express the conservation of energy and momentum us- 
ing T. Consider a cubical fluid element with each side of length /. Then amount of 
energy-momentum coming in must equal the amount of energy-momentum going out 
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of the cube. That is, 

(2.46) A /3Ta ° = ( Ta/ V = 0) - = 0} , 

/3=1 

where l 2 T aa [x u = 0) is the rate of flow in the x u direction across the surface of 
constant x a at x u = 0, and —l 2 T acr (x cr = I) the rate at x a — I with the minus sign 
indicating that the flow is coming in from the opposite direction of that at x a = 0. 
Now we divide each side of equation (|2.46|) by I 3 , and take the limit as I — > to find 

(2.47) T aP p = 0, 
where we have used the definition of the derivative 

(2.48) ton ^ = ")-^ = = Bp*, 
i^o I Ox 13 

In our Lorentzian frame, 

(2.49) = = 0, 
and since this is a tensor equation we can say 

(2.50) VpT a/3 = 

is the covariant expression for conservation of energy and momentum. 

2.3. The Einstein Equations 

The idea behind the Einstein equations is that the sources of the gravitational 
field determine the metric [|12|| . Classically, this is given by Poisson's equation 



(2.51) -A$ = AttQp, 

whose solution is given 

*(x) = I ^iy, 

V l x ~ y| 

where Q is Newton's gravitational constant. The function $ is the Newtonian gravi- 
tational source due the mass density p. This leads to the question, what quantity acts 
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as the source for the gravitational field? Since the stress-energy tensor T contains 
the mass-energy density p, the component T 00 , and is covariant we postulate T as 
the source for the gravitational field. Then to derive the relativistic generalization 
of equation ( |2.51| ), T should be coupled to a symmetric (0,2) tensor whose covari- 
ant derivative vanishes, and created from the geometry of spacetime. Finally, the 
coupling should reduce to the Newtonian equivalent given in equation (|2.51|) when 
velocities and the gravitational field are sufficiently small. The Einstein tensor G was 
specifically designed to satisfy these criteria. Therefore, 

(2.52) G = kT, 

with k = 8ttQ/c a , which reduces to equation ( |2.51| ) in the Newtonian limit of low 
velocities, and weak gravitational fields. For a perfect fluid the Einstein equations 
are written as 

(2.53) G aP = R af3 - ]^g aP R = 8vr£ (pg a/3 + (p + p)u a u p ) = %nQT a(i . 

2.4. The Second Fundamental Form 

The shock- wave solutions of the Einstein equations considered in chapters [3] and ^ 
are hypersurfaces of the spacetime manifold. A hypersurface is a n — 1 dimensional 
submanifold of an n dimensional manifold. The second fundamental form describes 
how the hypersurface is embedded in the spacetime manifold M by recording how the 
tangent spaces of the hypersurface change over the hypersurface. The first fundamen- 
tal form is the metric tensor g. In this section we will derive the second fundamental 
form. 

Let £ be a hypersurface of M, and denote the induced metric on S by g, which 
is assumed to be non-degenerate. For every point p in S, there exist slice coordinates 
x = (x l , . . . , x n ) on a neighborhood U of p in M such that U D E is given by x n = 0, 
and (x 1 , . . . , x n_1 ) form local coordinates for S |8|]. 
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Since the second fundamental form is a measure of how a tangent space changes on 
a hypersurface as is moves from point to point, it should involve the the connection 
V on the surface. Therefore, we need the idea of an induced connection on the 
hypersurface, and how to describe vector fields on the surface. 

Given vector fields X and Y on M, we already have already have a way to measure 
the vector rate of change of Y in the X p direction by the vector field on M. 

We should also note that the connection V is unique, since one can show that the 
condition, equation fl2.15|) , V a ^ 7 = implies that V is unique |18| , theorem 3.1.1]. 



We want to somehow induce the differential structure on to our hypersurface. 

Any smooth vector vector field X on a submanifold always has a local extension 



into its ambient manifold ||19|| . In our case, where E is a hypersurface of M, this 
means given a smooth vector field X on E, for each point p e E, there exists a 
neighborhood U C E of p and a neighborhood V C M of p such that U C V, and 
there exists a smooth vector field X on V such that 

X\u = X\u. 

For each p in E the metric g splits T P M into the direct sum 

T p M = T p E © (TpE)" 1 . 



Therefore, if X e T P M, then 



where X* a " e T P E and X" or e (T P E) ± . 

TTie Induced Connection. Suppose E is a hypersurface of the spacetime manifold 
M, and X, Y are vector fields on E. If X, Y are local extensions of X, Y to M, then 
we define 



(2.54) V ji Y = (V x Y) 

where V denotes the connection on E. 



tan. 
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Definition 2.4. If X, Y are local vector fields on E, then 
II(X, Y) = V X Y - V X Y = (V X Y) ± , 



is said to be the second fundamental form tensor [10]. 



The tensor II is a local vector field on M normal to X 

Definition 2.5. Let n be a unit normal vector field on a hypersurface E of M. 
The (1, 1) tensor field K on S such that 

(x(x),r) = (//(x,r),n) 

for all vector fields I,y on S is called the second fundamental form. K is also 
sometimes called the shape operator | |10| ] . 



Here, K determines a linear operator K : T p S — ► T p S at each point pGS [[10 

Lemma 2.6. Tne second fundamental form, 
(2.55) if = -V x n 

at eac/i point pGS, and if : T P S — > T P S. 

Proof. Since (n, n) is constant, ( , ) satisfies the Leibnitz rule, and is symmetric 
we have that (V^n, n) = for X G T P E. Thus, V^n is tangent to E for all X e T P S. 
Now, for all vector fields 7 on E, 

(K(X),Y) = (II(X,Y),n) = (V x y-(V x Y) ta ,n) = (V x Y,n), 

and since 

(V X Y, n) + ( Vxn, Y) = V x (Y, n) = 0., 

we have that 

(K(X),Y) = -(V x ii,Y). 

□ 



CHAPTER 3 



Lightlike Shock- Wave Solutions of the Einstein Equations 

The central results of this dissertation are contained in this chapter. The main 



result, Theorem |3.9| , is a set of equivalent conditions that imply conservation of energy 
across a surface in which the metric is only Lipschitz continuous. The other pertinent 
result, Theorem [3.24| , is a set of conditions which are equivalent to conservation of 
energy across a spherically symmetric hypersurface. In the case when the hypersurface 



under consideration is non-null, both Theorems 3.9 and 3.24 were first proved by 



Smoller and Temple in [|14| . The same results can also be found in |4[ , and |[L6| . Our 
contribution here will be to generalize Smoller and Temple's results to include the 
lightlike (or null) case. 

There are two main difficulties in the lightlike case. The first is that the met- 
ric restricted to the lightlike surface is degenerate. We will deal with this issue by 
considering the problem in the context of the whole spacetime manifold where the 
metric is not degenerate. The other, more problematic difficulty, is that the second 
fundamental form, defined by K = -V^n where X is a vector tangent to the surface 
and n is the normal to the surface, cannot be used to describe the dynamics of the 
surface in the ambient spacetime. This is because K measures the change in the nor- 
mal vector n in the direction tangent to the surface, but in a lightlike surface n lies in 
the tangent space of the surface; thus K can no longer give geometrical information 
about how the surface relates to the ambient manifold. We rectify this, using an idea 
of Barrabes and Israel given in [jlj , by defining a generalized second fundamental form 
/C = — VxN where N is a vector transverse to the surface. 

In the first section we define the generalized the second fundamental form /C, and 
construct a coordinate system analogous to Gaussian normal coordinates so that in 
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this coordinate system the components of /C can be written in terms of the derivative 
of the corresponding metric components in the direction of N. In the second section 
we state the main theorem of this dissertation. Then in the following section we give 
the supporting lemmas and theorems which will yield the proof of the main theorem. 
We finish off the chapter by stating and proving the theorem involving matching 
spherically symmetric metrics across a hypersurface. 

3.1. Generalizing the Second Fundamental Form 

Let M denote a manifold equipped with a metric g with fixed Lorentzian signature 
( — !-•••+), and E denote a hypersurface, which may or may not be null, that divides 
M into two regions M L and M R . Locally, we define E by ip(y) = 0, where ip is a 
smooth function such that 

(3.1) n t dy l = ^dy^0 

for any coordinate system. Let g L and g R denote the metrics on M L and M R respec- 
tively, and the metric on M is given by g = g L U g R . Assume g L , and g R are smooth, 
that is at least C 2 , on each side of E with their derivatives uniformly bounded across 
E. 

The failure of the Second Fundamental Form in the Lightlike Case. 

Given a point, p in E let T P E denote the tangent space of E at p. A vector X is in 
T P E if 

(3.2) (n,X)=g a pn a XP = 0. 

If E is a lightlike hypersurface, then (n, n) = 0, which implies that n is also in T p E. An 
unfortunate consequence of this is that [K] = is always true, where [ • ] = /C L — K, R 
denotes the jump in /C across E, and hence does not yield any information about how 
E is embedded in M. Therefore, we cannot use the standard second fundamental 
form to determine the necessary and sufficient conditions for a null hypersurface to 
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be a shock surface or a surface layer. This failure of the second fundamental form to 
yield conditions on the existence of lightlike shocks is a consequence of n not being 
transverse to E; therefore any alternative to K, designed to extract this information, 
should involve a vector transverse to E. 

Generalized Second Fundamental Form. We will generalize the second fun- 
damental form given by equation (|2.55|) so that it will still be defined in the case of 
lightlike hypersurfaces. We will use an idea introduced by Barrabes and Israel in [|TJ. 
Even though they use a scalar version of the second fundamental form and we are 
not, the idea is still the same. 

We begin by choosing a vector N transverse to E, that is, a vector not in T P E. 
Then we replace the second fundamental form K(X) = -Vjn with what we will call 
the generalized second fundamental form 

(3.3) K{X) = -V X N, 

where X is in T p E. If we compare this with the definition of the second fundamental 
form in equation ( |2.55| ), then we notice that /C depends on the vector N that we 
choose just as K depends on the normal vector n. To make K, well-defined we need 
to place certain restrictions on which transverse vectors N can be chosen to define /C. 
We need N to be continuous across E, so we require for all points p in E 

(3.4) [(N,X a >] = [iV Q ] = 0, 

where {X a }"~J is a basis of T P E. Furthermore, the length of N on each side of E must 
coincide; hence 

(3.5) [(N,N)] = 0. 

The transverse vector chosen under the conditions ( |3.4| ) and ( |3.5| ) is not unique. Any 
vector N satisfying equations Q3.4Q and ( |3.5| ) is invariant under the transformation 



(3.6) N^N' = N + A a X a , 
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where, again, {X a }"~J is a basis of T p E, and A a is an arbitrary function. Thus, we 
are free to choose N so that 

(3.7) (N,n>=r^0, 

everywhere on E, where rj is any non-zero constant we like. Now, under the transfor- 
mation (|3.6| ), JC transforms as 

(3.8) JC b a X a — ► JC b a X a - X c T b ac X a . 
Equation (|3.8|) comes from the computation 

- (Vx (N + A%)) 6 = JC b X a - (X\ a + X c T b ac )X a = JC b a X a - \ c T b ac X\ 



Equation ( |3.8| ) tells us that JC does not transform as a tensor, since \ k T J ik does not. 
Since we are seeking a tool to describe how E in embedded in ambient spacetime 
that is coordinate independent, this presents a slight difficulty This is rectified by 
the insight of Barrabes and Israel in |]J that [JC] is a tensor. In the lemma below 
we will show that [JC] is invariant under the transformation (|3.6| ), thus [JC] does not 
depend on the particular transverse vector N that is chosen other than that it must 
satisfy the the jump conditions (|3.4| ) and ( |3.5| ). The reason for the invariance is that 
the metric g and its tangential derivatives are always continuous across S, but its 
transversal derivatives may not be. 

This means we can now say that JC is a true generalized form of the standard 
second fundamental form K when we compare their jumps across E. That is, if E is 
non- light like, then we choose n, or another suitable vector, as our transverse vector 
we get [K] = [JC], and in the limit, as E goes lightlike, [JC] remains well-defined; hence 
[JC] is the tensor we seek to determine how E is embedded in spacetime even when E 
is light like. 

Lemma 3.1. The jump [JC] = (JC R — JC R ) across E is independent of the choice of 
N, and is a tensor. 
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Proof. Lowering the index j in equation fl3.8| ) we get 
(3.9) KL ab X a — > JC ab X a - \ c g dh Y d ac X a . 



Since 



9db^ac — ^(—9ac,b + 9ba,c + 9ca,b) 



only involves derivatives in the tangential direction, and 

(3.10) [)C ab - \ c g lb T l ac ]X a = [IC ab ]X a - A c [g lb T l ac ]X a , 



under the transformation (3J3), we have 

(3.11) [JC ab ]X a — > [JC ab ]X a - A c [g db T d ac ]X a = [JC ab }X c 



□ 



Modified Gaussian Skew (MGS) Coordinates. We construct a new set of 
coordinates which will make computation of the Riemann and Einstein curvature 
tensors more manageable. 

One important ingredient in characterizing solutions of the Einstein equations 
across an interface, first formulated by Israel H, and also used in Smoller and Tem- 



ple's work ||14|, 15, [16], is that in Gaussian normal coordinates you can write the 



second fundamental form in terms of the metric as 

1 

K a b — —-gabfl, 

where represents the 0th coordinate, whose corresponding basis vector is the normal 
to E. In Gaussian normal coordinates, values of [K] can be used to easily compute 
the corresponding jumps in the Riemann, Ricci, and Einstein tensors. In the lightlike 
case, however, Gaussian normal coordinates are undefined since the normal vector 
n is a tangent vector. To rectify this we employ Modified Gaussian Skew (MGS) 
coordinates. 

A MGS coordinate system is constructed in an manner analogous to Gaussian 
normal coordinates by replacing the vector normal to £ with a transverse vector that 
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satisfy the jump conditions (|3.4j ) and (|3.5| ). The property that the Oth coordinate 
vector is orthogonal to the other n — 1 coordinate vectors is lost in a MGS coordinate 
system, but such a coordinate system can be defined even when E is a null hypersur- 
face. However, to deal with the degeneracy of the vector subspace T p E we will place 
certain restrictions on the basis of T P E we choose. These additional restrictions will 
define MGS coordinates, and in doing so, the metric, in MGS coordinates, will be 
almost diagonal. 

We model the construction of MGS coordinates after the construction of Gaussian 



normal coordinates by Smoller and Temple from ||14| , |16|| . One can construct MGS 
coordinate as follows: 

• We first choose coordinates in the surface. Let T P E denote the tangent space of 
E at point p G E. The surface coordinates y 1 , . . . , y n ~ 1 will be chosen to coincide 
with a particular set of basis vectors of T P E. First, choose y 1 so that ^p- = n. 
The remaining y\ i = 2, . . . , n — 1, are chosen so that = {^-,-^1^ — for 
i 7^ j. We will show below that ga > 0. Furthermore, we will also place the 
restriction on each y l that 

(3-12) (^' N H' 
for i = 2, . . . , n — 1. 

Notice that we can always choose such a coordinate system. First, n is in 
T p E, and is non-zero, so we can choose an orthogonal basis incorporating n. 
Also recall that once N is chosen we are free to modify it via a transformation 
of the form given in equation ( |3.6|) . Thus we can always transform 

(-3) N^N-g(£,N)£ 



in order to satisfy equation (|3.12j ). 



For each p 6 E, let 7 P (s) denote the geodesic which satisfies 

7 P (0) = p, 7 P (0) = N, 



3.1. GENERALIZING THE SECOND FUNDAMENTAL FORM 37 

where N is a vector transverse to E satisfying the jump conditions ( |3.4| ) and ( |3.5| ). 

• Assume N points into the right hand side of E, for convenience. 

• Choose coordinate w° so that if j p (s) = q, then w°(q) = s. As a consequence, 
w° < on the left hand side of E, and w° > on the right hand side of E. 

• The coordinates w 1 , . . . , w n ~ l are given by w l (p) = y l {p) of p G E, and w l (q) = 
w l (p) if and only if q = 7 p (s), for some p and s, where i = 1, . . . , n — 1. 

By construction, we can justify the following lemma. 

Lemma 3.2. In MGS coordinates, 
(3.14) ds 2 = g 00 d{w ) 2 + 2r] dw^dw 1 + g ii dw i dw\ 

where i = 2, . . . , n — 1, and rj = g(N, n). 

We now show how an MGS coordinate system relates at other coordinate system. 

Lemma 3.3. There exists a C 1,1 transformation from any coordinate system to an 
MGS coordinate system. 

Proof. The hypersurface E is given by ip{y) = in y-coordinates. Using a 
smooth transformation we can use slice coordinates to define E by ?/ = 0. We 
assume that E is lightlike, hence we can write the normal to E in terms of vectors in 
T p E for any p in E, so that 

n = na ^ a ^°- 
Let n P {t) denote the geodesic which satisfies 

M°) = P> Ap(°) = n - 

Then we can smoothly transform the slice coordinates y to a new set of slice coor- 
dinates u, where u 1 = t and u a = y a for a ^ 1, so that E is given by u° = 0. In 



^-coordinates we are free to use the transformation (3.13) to adjust so that 
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but still have (N, n) = rj. Then our MGS coordinates can be written as w = 
(w°, . . . , w n ~ 1 ) = (s, t,u 2 , . . . , u n_1 ). Now for q not in S, but on the geodesic 7 P (s) 
and in the coordinate neighborhood, we have that w(q) = (s, 0, u l {p), . . . , M n_1 (p)) is 
a smooth function on each side of S. Therefore, it only remain to check the continuity 
of the derivatives across at s = 0. Indeed, 

dw» b ' 



and 

du a 



dw° 

where N a denotes the components of N in ^-coordinates. This comes from the vector 
transformation law: 

N = _d_ = dv^_d_ = NO , d 



dw° dw° du a du a ' 
Since N is continuous across £ it follows that the derivatives are as well. □ 



Remark 3.4. Notice that it is possible for N to be null, which would mean goo — 
in MGS coordinates. We will see in the proof of the lemma below that the span of 
the vectors N and n forms a two dimensional Lorentz vector space. Since any two 
dimensional or greater Lorentz vector space contains two linearly independent null 
vectors, it is possible for both N and n to be null. For more information regarding 
Lorentz vector spaces see [jl^, Chapter 5]. 



Lemma 3.5. In MGS coordinates, the metric components gu > fori = 2, . . . , n— 
1, and hence the coordinate basis vectors can be normalized so that gu = 1 for i = 
2,... ,n-l. 

Proof. Let T p M denote the n-dimensional tangent space of the spacetime man- 
ifold Mat a point p, and let peS. Now choose MGS coordinates on a neighborhood 
U containing p. Then the coordinate basis of T P M is (N, n, e2, ... ,e n _i), where 
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e 2 , ■ ■ ■ , e n _i are orthonormal with respect to each other, and are orthogonal to both 
N and n. Thus, in MGS coordinates at p the metric g is given by the matrix 



(3.15) 



1 T] 
7] 

g 22 












\ 







\ u 

where rj is a constant, and ga — ± 1 for i = 2, . . . , n — 1. Notice that for a transverse 
vector satisfying the jump conditions (|3.4j ) and Q3.5|), it is possible for g(N,N) = 
0, ±1, which would make goo = 0, ±1 in MGS coordinates. However, we can always 
"adjust" N by a transformation of the form ( |3.6| ) to make goo — 1- If <?(N, N) = 0, 
then let N -> N + n/2rj, and if g(N, N) = -1, then let N -> N + 11/77. Therefore, 
we can say g 00 = 1 in any MGS coordinate system. 

Now, since T p S is a degenerate subspace of the Lorentzian vector space T p M, 
and gn = g(n, n) = 0, we must have ga > because a degenerate subspace of a 
Lorentzian vector space cannot contain any timelike vectors, and also cannot contain 
more than one linearly independent lightlike vectorQ. □ 



Remark 3.6. Using the idea from the above proof we can construct a coordinate 
system so that the metric g is of the form diag(— 1, 1, . . . , 1), which we will use later 
in this section as a candidate for an inertial coordinate frame. We construct these 
coordinates, denoted by u 



u°, . . . ,u n 1 ), on the neighborhood U of p in S from 



MGS coordinates w = (w°, . . . , w n 1 ) on U of p. From Lemma p.5| we can choose 



w l so that g{w l ,w l ) = 1 for i = 2, ... ,n — 1. Now we let u 1 



w ,u 



w l for 



i = 2, . . . , n — 1, and choose u so that 







1 d d 



du° \rj dw 1 dw° 



(n — i]~N)/i]. 



^ee [1C, Chapter 5] for more information on vector subspaces of Lorentzian vector spaces. 
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Then, in u-coordinates, g 00 = g ((n — r]N)/r], (n — r]N)/r]) = — 1, g n = g(N, N) = 1, 
and go\ = g ((n — r]~N)/r], N) = 0. Thus, in w-coordinates, the metric g has the form 
diag(-l,l,...,l). 

The Main Result of the Generalized Second Fundamental Form. The 

following lemma shows that the components of the generalized second fundamental 
form can be written in terms of the derivative of the corresponding component of 
the metric g in the N direction. As a consequence, from values of the jump in the 
transverse second fundamental form, [K], we will be able to determine the values of 
the jumps of the Riemann curvature tensor [i?£Lg], and hence the Ricci tensor [R a p], 
the Ricci scalar [R], and Einstein curvature tensor [C7 Q( g]. 

Lemma 3.7. In MGS coordinates 
(3.16) IC ab = --0a&,o, 

where a, b = 1, . . . , n — 1. 

Proof. For any vector field X 13 on the surface S we have 



where j3, r, a = 0, . . . ,n—l. The last equality comes from the property that in MGS 
coordinates the components of the transverse vector N are N a — 0, a — 1, . . . , n — 1, 
and iV° = 1. 

Since g a o is constant for all a = 0, . . . , n— 1, the Christoffel symbol in equation 



(3.17) 



= (VxN)^A^ 

= (fiTp + I^) X? 




( p. 17] ) can be written as 



(3.18) 



/30 — 



2 id™ {—9f30,a + 9a/3,0 + dOa^}} ~ -g^dabfl- 
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where a, b — 1, . . . ,n — 1. Combining equations ( |3.17| ) and Q3.18D gives 




Lowering the index r in equation (|3.19| ), where 



(3.20) 



we get the desired result 



(3.21) 



fca.b — ~7:9ab,0- 



□ 



This result will be the key tool we will use to prove the Main theorem of this chapter, 
and of this dissertation! We end this section with a corollary of the above lemma 
which will we use to prove part of our main theorem. 

Corollary 3.8. The metric components of g = g L Ug R in MGS coordinates are 
C 1 functions of the coordinate variables if and only if [/C] = {JC R — /C L ) = at each 
point on the surface S. 

Proof. The jumps in derivatives of metric components in the direction of the 
surface £ are always zero, that is, g a f3, c for a, /? = 0, . . . , n — 1 and c = 1, . . . , n — 1, 
since [g a p] = on E. Also, g p is constant for each (3 = 0, ... ,n — 1 Then by 
Lemma [T7|, [K. a b] = —^[g a bfl] = implies that g is C 1 for each coordinate. □ 



The chief result of this dissertation is the theorem stated below which classifies 
all Lipschitz continuous solutions of the Einstein equations across a hypersurface. 
Smoller and Temple first proved this result in the case of non-lightlike hypersurfaces 



the case of a lightlike hypersurface using the generalized second fundamental form, 
developed in the previous section. 



3.2. The Main Theorem 



in [14], and can also be found in |3) and [16]. We extend their results here to include 
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Theorem 3.9. Let E be a hypersurface in spacetime, and let JC denote the gen- 
eralized second fundamental form on S. Assume that the components g a @ of the 
gravitational metric g are smooth on either side o/E, (continuous up to the boundary 
on either side separately), and Lipschitz continuous across E in some fixed coordinate 
system. Then the following statements are equivalent: 

(i) [K] = at each point o/E. 

(ii) The components of the curvature tensors R^p o,nd G a p, viewed as second order 
operators on the metric components g a p, produce no delta functions on E. 

(Hi) For each point p in E there exists a C 1,1 coordinate transformation defined in a 
neighborhood of p, such that, in the new coordinates, (which can be taken to be the 
MGS coordinates for the surface), the metric components are C 1,1 functions of these 
coordinates. 

(iv) For each p in E, there exists a coordinate frame that is locally Lorentzian at p, 
and can be reached within the class of C 1 ' 1 coordinate transformations. 

Moreover, if any one of these equivalencies hold and E is non-null, then the 
Rankine-Hugoniot jump conditions, [G]pN^ = 0, which express the weak form of 
conservation of energy and momentum across E when G = kT , hold at each point on 
E. 

Lastly, if we add the condition that the second derivative of the the metric inner 
product on T P E + is continuous, where T P E + denotes the spacelike subspace of T P E ; 
then the Rankine-Hugoniot jump conditions, [G^N 13 = 0, also hold for a null-surface. 
In MGS coordinates, the condition that the second derivative of the the metric inner 
product on T P E + is continuous is equivalent to the condition [gu,oo] — 0, for t = 
2,... 

3.3. Supporting Results 

The proof of this theorem involves series of lemmas and theorems which are de- 
veloped below. 
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The Rankine-Hugoniot Jump conditions. The objective of this section is to 
write down the Rankine-Hugoniot jump conditions in terms of our generalized second 
fundamental form, which will give the criteria for when the jump conditions hold. 
These jump conditions, 

(3.22) [0\pN p = 0, 

express conservation of energy and momentum on the shock surface £ when G = kT 



in the weak sense W%> We begin with a lemma enables us to compute the 



connection coefficients in terms of the generalized second fundamental form. 

Lemma 3.10. In MGS coordinates the components of the connection coefficients, 
for a metric g at a point p m £, can be written as 

(3.23) r ab = r ab + ^X, 

(3-24) Tl = -JCl, 



where 7 = 0, . . . , n — 1„ a,b,s = 1, . . . , n — 1, and Tl b indicates the summation of 
the indices on the metric components goes from 1 to n — 1. 

PROOF. Writing the connection coefficients in terms of the metric components we 
see that 



(3.25) 



~g JS {-g a b, s + g sa ,b + gbs.a} + ^g l0 {-gab,o + goa,b + gm,a} 
Kb - \g l0 g a bfi 



g^K 



ab ■ 
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which is equation ( |3.23| ). We use the same method to get equation ( |3.24| ), that is, 



Flo — 2# 7fT {~ 9a0,a + gaafi + gOa,a} 
(S- 26 ) = ^g JS {-9aO,s + 9sa,0 + 90s,a} + ^fi ,7 °{-fi'aO,0 + 90a,0 + 900,a} 

= ^ S 9sa,0 = -JQ. 



□ 



With the use of the above lemma we can now write the components of the Riemann 
curvature tensor in terms of the components of the generalized second fundamental 
form. 

Lemma 3.11. In MSG coordinates the components of the Riemann curvature ten- 
sor can be written as 

(3.27) Rl jb = {JC ab;j - K aj , b ) + {V h K, aj - K]K ab ) 



(3.28) 



and 



(3-29) Rl b = ^l b -n,r 



Recall that the Greek indices go from 0, 



, n — 1 and the Latin indices go from 
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Proof. For equation ( 3.27|) we will first consider three equations that will be 



pertinent to the computation. Using Lemma 3.10, we can write 



(3.30) 



— g^° {fc-ab,j — K-vbFaj ~ (fcaj,b ~ ^aj^ab)} 

— g^ (fcab;j — K-aj;b) , 



(3.31) 



and 



(3.32) 



p7 -ps _ p7 | -ps 
1 sb L aj — L sb I 1 aj 



f 7 v s — r 7 v s _i_ n s0 r 7 jt 

1 sj L ab ~ L sj 1 ab ' g L sj'^abi 



p7 pO 
1 0j L ab 



sy 



-q(r° ab + g 00 K: 



ab 



1 ab^-j ~ ^ab^j 



p7 pO _ _p0 j^7 _ ^ u~y 
1 Ob 1 aj — 1 aj'H '^aj'^b • 



Using Lemma |3.10| , and equations ( |3.30| ), ( |3.31| ), ( |3.32| ) we have 



c>7 — F 7 _ F 7 _i_ F 7 F CT _ F 7 F CT 

-^ajb ~ 1 ab,j 1 aj,b ' 1 aj L ab 1 ab 1 aj 

— F 7 _ F 7 _i_ F 7 F s _ F 7 V s _L F 7 F° _ F 7 F° 

~~ 1 o6,j 1 aj,b 1 sj 1 ab 1 sb 1 aj ' 1 Oj 1 ab 1 Ob 1 aj 



(3.33) 



Kb, + g" ^ - r;. b - g ^K aj , b - r;. (v] b + <f X 6 
+ ^Ib frjj + g lQ K-sj^ + r^r^ — r^rj^ 

g' 70 {fcabj - )C sb r S aj - (JC a j,b - KsjKb)} 

, p7 f s _ F 7 V s _i_ F 7 _ F 7 _l_ F 7 F° _ F 7 F° 
1 sj 1 ab 1 sb 1 aj ' 1 ab,j 1 aj,b ' 1 Oj 1 ab 1 Oft 1 aj 

g ^ (iCa^ - K aj , h ) + (jqKaj - qic ab ) + (i%iq - f °jq 



+ g s0 M^ab-riJc a A+R: jb , 



whprp R 7 — F 7 — F 7 -I- F 7 F s — F 7 F s 
wneie _n, ajfe — i ab - i aj . )6 -t- 1 sj L ab i aj . 
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For equations (|3.28|) and (|3.29|), we also use Lemma |3.10| to write 



E>7 



(3.34) 



and 



r 7 

1 aO,j 



-K? ■ 



p7 I p7 -ps _ p7 -ps | p7 -pO 
1 aj,0 1 sj 1 aO 1 sO 1 aj "r 1 Oj 1 aO 



p7 pO 
1 00 1 aj 



r 7 • - <7 7 X;,o - >q r 7 + g ^K 



+ K 1 ( Kj + 9 S °£aj ) + 



r 7 



p7 I p7 per p7 per 



Oj,b 



ah 1 - Oj 



(3.35) 



r 7 - k 1 



□ 



Remark 3.12. For any quantities A, and B such that both [A] = and [5] = 0, 
the jump in the product of the quantities is 

[AB] = A L B L - A L B R + A L B R - A R B R = A L [B] + B R [A) = 0. 

The next three lemmas are more of the technical variety which we also need to 
facilitate computations to be made later on. 



LEMMA 3.13. If [JC a b] = for every point in S, then [JC a b-j] = for every point in 



E. 



(3.36) 



Proof. Using Lemma |3.10| we have 

Since [K, a b] — for every point in £ implies that [/C a b,j] = for j 
[r^-] = 0, the result holds by equation ( |3.3tj| ). 



fcabTlj 



. , n — 1, and 



□ 



Lemma 3.14. // [K] = 0, then [R 1 *] = 0. 
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PROOF. If [K\ = 0, then by lemmas |3.11| and p,13| , remark |3.12| , [r 
[Rlj b ] = 0, we can conclude both 

(3.37) K,J=0 and [fljj = 0. 

Therefore, if [/C] =0, then 

(3-38) = <T K i6 ] = <T KJ + 9 a ° [RU = 0. 



7 1 
ab\ 



0, and 



□ 



Lemma 3.15. If [K] = 0, then in MGS coordinates 

[K\] = 9 a °[9u, t 

for a = 1, . . . , n — 1, and t > 1. 



;,00 



Proof. First, 



(3.39) 



pat „tcr pa 



-g u R a 



tto> 



since for t > 1, g u are the only non-zero components. We now compute R\ m in MGS 
coordinates by using lemma |3.10| which yields 



pa 



■pa I "pa ps pa ps _j_ pa pO pa pO 

1 tt,0 st £0 sO 1 ii "T 1 0t L tO ~~ 1 00 1 



to,t 



(3.40) 



When [/C] = we have 



-St 



+ lCUT s tt + g s0 !C tt )+}C?}C°. 



pc- j = <r [/c rt ,j = o, 



pa 
1 tt,0 



[fl'tt.o] a = ' 



and the remaining terms are zero, when [K] =0, except 
(3.41) [ 5 a0 /C« >0 ] = g°°[g m 

when 



:,00J 



□ 
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The next result gives us a nice way to express the components of the Einstein 
curvature tensor in terms of the components of Riemann curvature tensor. We repro- 



duce the statement of the following lemma, and its proof, given in [f|], ||14|| , and |T6 
These same computations can also be found in [^, section 14.2]. 

LEMMA 3.16. The components of the Einstein curvature tensor can be written as 

(3.42) GS = -E i2 IH« « = V--,n 
and 

(3.43) G% = R \7t\> a ^P- 

The indices inside \ ■ \ are always taken as an increasing sequence. 

Proof. We begin by raising the index of the components of the Einstein tensor 
defined in section |2.1| , and given by 

from which we find 

(3.44) Gp = Rp — \g$R = B$ - ±fi$R. 
Since R a J v is antisymmetric in err and also in fxu, it follows that 
(3-45) R a a = RZ = J2 R ^ = R S- 

Notice that we do not sum over the index a since it represents the component, and 
not a summation index. The Ricci Scalar can be written as 

(3.46) R = RZ = 2R\Z\. 

Combining equations ( |3.45| ) and ( |3.46| ) we have equation ( |3.42| ), 

not _ o pl a H pl CTT l _ pl CTT l 

— Ztt, i — it, i — — > it, i. 

a \ al/ \ \°" r \ / j \ar\ 
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Furthermore, we also find equation ( |3.43| ) by the computation 
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^(3 — n f3 — n f3r 



Epkl 



□ 



Now we are finally able to directly prove a part of our main result, Theorem p79 



The following theorem gives the criteria for when the Rankine-Hugoniot jump condi- 
tions hold across S. 

Theorem 3.17. If both [g tt ,oo] = for t > 1 in MGS coordinates, and [K] = 
at each point ofH, then [G a /3}N a = 0. The condition [gtt,oo] — for t > 1 in MGS 
coordinates means the metric is C 2 in coordinates not associated with the transverse 
vector N and the n to the surface. 



Proof. In MGS coordinates N° = 1 and N a = for a = 1, . . . n — 1, therefore 



(3.47) 



When (3 = 0, we have 



(3.48) 



[G ap }N a =[G op }=9o [G ,}+g i[Gl} 



[Go/3] — [Gqo] — 9oo [Gq] + <?oi [G ] 



(3=1, 
(3.49) 



[Gop] = [G 01 ] = 9i [G° ] + g n [Gl] = g 10 [G° ], 



and when (3 = t > 1 we have 



(3.50) 



[Gop] — [Got] — gio[G ]. 



See equation (|3.15 ) for the components of g in MGS coordinates. By lemma |3.16 
have 

1 4£[« 



we 
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and, by lemma |3.14| , [JC] = implies that [R^] — for all s, t — 1, . . . , n — 1, hence 



[G° ] = 0. Also, by lemma \TT§ 

[G a ]- 



E <' 



|at| 



Now, lemma |3.15| , [/C] = implies that 



ra = <? a v,oo], 

for all t = 2, . . . , n — 1. Then [gt^oo] = for £ > 1 implies that [Gg] = 0. Therefore, in 
MGS coordinates, 

[G aP ] N a = 0. 

□ 

Curvature Tensor Components as Delta Functions. The physical condition 
for the hypersurface S to be considered a shock surface, when G = kT, is that there 
exist no delta function singularities at a point p in E in the components of the Einstein 
curvature tensor G. 

What does this mean? When considering weak solutions of the Einstein equa- 
tions G = kT across a discontinuity if there exists a delta function singularity in a 
component of the Einstein tensor, then this singularity manifests itself in the coupled 
stress energy tensor T as an infinite spike in the mass/energy density of matter at 
the discontinuity || . In this case the hypersurface of discontinuity is referred to as a 
surface layer. On the other hand, if there exist no delta function singularities in the 
Einstein tensor, then there is merely a discontinuity in the mass/energy density, and 
the surface is considered a shock. 

We will show below that there exist no delta function singularities in the Riemann 
curvature tensor if and only if [K] = for every point on E, and hence will give a 
necessary and sufficient condition for a hypersurface to be considered a shock surface. 

To see this, notice that each component of the Einstein tensor can be written 
in terms of the metric components g a p, its first and second weak derivatives, and 
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the components of the inverse of the metric g a/3 . Since g a p and g a ^ are Lipschitz 
continuous across £ it follows that the only possible delta function singularities con- 
tained in the the Einstein tensor are the second order derivatives, g a p^ri of the metric 
components. Also, notice that derivatives in the direction tangent to £ always can- 
cel out because [g a p] — at all points of E. This means that, in MGS coordinates, 
[9a/3,s] = for s = 1, . . . , n — 1, and the only possible delta functions are of the form 
9ab,oo- Therefore in MGS coordinates, when [K ab ] = \g a b,o) — 0, the metric g is C 1 
just as we saw in corollary |3.8| , and g a b,oo is n °t a delta function. However, when 
[K-ab] = [9ab,o] 7^ 0, then g abfl0 is a delta function. 

In the following lemmas we will show that components of the Riemann curva- 
ture tensor, and hence the components of the Ricci tensor, the Ricci scalar, and the 
Einstein tensor contain no delta functions if and only if [K] = 0. 

Lemma 3.18. In a MGS coordinate frame, with respect to the surface S ; the com- 
ponents of the Ricci tensor can be written as follows: 

(3.51) i?oo = —-^g ab 9ab,oo + lower order 0-derivatives, 

(3.52) R a0 = —-g l0 g a i,oo + lower order 0-derivatives, 

(3.53) R a b = -g a b,oo + lower order 0-derivatives. 

Furthermore, the Ricci scalar can be written as 

(3.54) R = —g s0 g to g s t,oo + lower order 0-derivatives. 
Proof. From equations ( |3.27|) ( |3.28|) and (|3.29|) we have 



(3.55) Rlj b = lower order 0-derivatives, 

1 
2 



(3.56) Rlj = —-g^gajflo + lower order 0-derivatives, 



(3.57) R^ b = lower order 0-derivatives. 
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Thus, using Lemma |3.10| we can write equation ( |3.51| ), 



p pa _i_ pO pa 

-n-OO — -n-OaO -"-000 ~~ -"-OaO 

"pa "pa _j_ "pa per pa per 

~~ 1 00,a 1 0a,0 ~r 1 era 1 00 1 crO 1 0a 



(3.58) = K,o - r: roa 



^ ^as _ v~a y~s 

-^9 9as,00 ~ ^ s ^a 

— -g as g as ,oo + lower order 0-derivatives, 



where we have used Tq = 0. Also, from equation (|3.56|) we have 



(3.59) R a0 = R s as0 + R° a00 = R s as0 = ~g s0 g as ,oa + lower order 0-derivatives, 



and from Lemma |3.11| and equation (|3.34|) we have 



(3.60) R ab = R s asb + R° a0b = ~5W>,oo + lower order 0-derivatives. 

Using the results above we can write 

R = g a<3 R a p = g ab R a b + g a0 R a o + g ob Rob + g 00 Roo 
= g ab R a b + 2 g a0 R a o + g 00 Roo 

= ^g ab g a bflo + 2g a0 (-^g W g a t,oo^ - ~9 ab 9ab,oo + lower order 0-derivatives 
— —g a °g t0 gat,oo + lower order 0-derivatives. 

□ 

Corollary 3.19. Assume g = g L U g R is smooth on either side of E ; and Lip- 
schitz continuous across S. Then in MGS coordinates the jump condition [JC] = 
exists at a point p in S ; if and only if, the curvature tensors R^g, R a fi, G aj 3, and the 
Ricci scalar R, viewed as second order differential operators in the weak sense on the 
metric components g al 3, produce no delta function singularities at the point p in S. 
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Proof. By applying lemmas |3.11| and |3.18| , we see that, in MGS coordinates, the 



components R^g, R a p, G a p, and the Ricci scalar R can be written in the form 

A{x)g ab fiQ + lower order O-derivatives, 
where A(x) is some function of the coordinates which is as least Lipschitz continuous. 



From Lemma 3/7 we know, in MGS coordinates, [K] = at a point p in E, if and 
only if, [g a bfi] = for all a, b = 1, . . . , n — 1 at the point p in S. We also know g a b,oo 
is not a delta function, if and only if, [g a b,o] = 0. This completes the proof. □ 



Now our goal is to generalize corollary |3.19| for any coordinate system, not just 



for a MGS coordinate system. To do this we will need the following lemma which 
considers how the Riemann curvature tensor, defined in terms of second order weak 
derivatives of the metric components g a p, transforms from one set of coordinates 
to another by a C 1,1 transformation. Recall by a C 1,1 transformation we mean a 
function whose first derivatives are Lipschitz continuous. The lemma was originally 
proved by Temple and Smoller in ||14| , and can also be found in |3| and in | 16| . We 



will reproduce their result here for convenience. 



Lemma 3.20. Let 1Z = R^ = L[g] denote the components of the Riemann curva- 
ture tensor in x- coordinates where L is the second order linear operator on the metric 
components g a p which defines R^. Similarly, let 71 = R^s = L[g] denote the com- 
ponents in y- coordinates which are related to x-coordinates by a C 1,1 transformation. 

If H, is a weak solution of1Z = L[g] in x-coordinates, then 7£|^ is a weak solution 
ofTZ = L[g] for any coordinate system y related to x by a C 1 ' 1 transformation. Here 
we have used the short-hand notation 

dx A dx^ dx u dx^ dy a 
dy~ = ^dyPdyidytdx*' 

and multiplication by a function is taken in the weak sense. 
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Proof. Let g be smooth and ip be an arbitrary smooth test function with compact 
support. Furthermore, let 

/ L\g]<p = [ L*\g,<p], 
Jr 4 Jr 4 

where L*[g,ip\ is defined as the expression obtained from L[g] by integrating the 
second order derivatives in g once by parts. Now, in any coordinate system, L is 
given by 

(3.61) 

= L[g] = — £ + lower order O-derivatives 

+ lower order O-derivatives 

= 9 Xa {-9vt,av + 9ic,nv + g^M - 9™m) + lower order ^-derivatives. 
Thus, L*[g,(p] is composed of the metric components g^ u , the test function tp, their 
first derivatives, and the inverse metric entries g^ v '. Therefore, L*[g,ip\ is integrable 
over M 4 for any Lipschitz continuous metric and any Lipschitz continuous test function 
<p with compact support. 
Now, suppose 

{K,<p) = I -top = I L*\g,ip], 
jr 4 Jr 4 

for all Lipschitz continuous test functions p with compact support, that is, 1Z is a 
weak solution of 1Z — L[g\. Let 

_ _ _ dx^ dx u dx 

9 = 9*0 = 9^-^0^ = 9-^- 

If || is Lipschitz continuous, then L* [g(dx/dy), ip] is bounded for any Lipschitz con- 
tinuous test function. 

Suppose g is an arbitrary, non-degenerate, Lipschitz continuous metric, and ip an 
arbitrary Lipschitz continuous test function. Furthermore, suppose the coordinates x 
and y are related by a C 1 ' 1 transformation, that is, |^ e C ' 1 Let g e a p denote the 



3.3. SUPPORTING RESULTS 55 

approximation of g a p by a smooth function, and also let x e (y) denote the approxima- 
tion of the coordinate map x(y) where x e (y) is smooth, and has a smooth inverse^]. 
These approximations can be chosen so that 



9 a p 



g af3 in C 0,1 , x e (y) — ► x(y) in C 



1.1 



ffo)-gfo)mC- and ^ *(,) i n 



Consequently, 



Now, if we define 



g< = g<p- ^ g , and ^ — # in C ' 1 . 
ax 6 



(3.62) 



and 
(3.63) 

then, by definition, 
(3.64) 



^,,dx e 
dy 



TV = TV 



dy_ 

dx e 



L* 



f dx e 
dy 



for all test functions <p. From this equation we can say TV^- is the Riemann curvature 



tensor for g e ^- Now, because every function in equation (|3.64j) is sufficiently smooth, 
the equation holds, if and only if, 



(3.65) 



Since g e — > g as e — > in C 0,1 , if follows that 



/ L*[g e ,cp) — > / L*\g,<p] as e -> 0. 
</r 4 </r 4 



2 In our case here, it is sufficient that if u is a continuous function on a bounded set, then there 
exists a smooth function u e such that ||u e — u||oo — ► as e — > 0. For more on approximations by 



smooth functions see §, Section 5.3] or |l3| , page 62]. 
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Therefore, from equation ( |3.65| ), TZ e — > TZ as e — > in the weak sense. In the same 



way, lZ e — > 1Z, and thus, by equation (|3.63|) , we conclude 1Z = in the weak sense, 



which completes the proof. □ 
With Lemma |3.20|, we are now able to prove the result in corollary |3.19| for any 



coordinate system with our Lipschitz continuous metric g = g L U g R . 

Theorem 3.21. Assume g = g L U g R is smooth on either side o/E, and Lipschitz 
continuous across E. Then the jump condition [JC] = exists at a point p in E ; if 
and only if, the curvature tensors R^g, and G a p viewed as second order differential 
operators in the weak sense on the metric components g a/ 3, produce no delta function 
singularities at p in E. 

Proof. By corollary |3.19| , the theorem holds in MGS coordinates. Now, for any 
metric g = g L U g R which is smooth on either side of S, and Lipschitz continuous 
across E, a transformation from our MGS coordinate system to another coordinate 



system must be least C > by lemma |3.3| . Since the transformation is an invertible 



C 1 ' 1 function, by Lemma |3.20| , the theorem holds for any coordinate system. □ 

Relating to a C 1 ' 1 metric by a C 1 ' 1 transformation. Here we prove a result 
that relates values of [JC] to the smoothness of the metric. In particular, the condition 
that [JC] = for each point p in E is equivalent to the existence of a coordinate 
system at each point in E such that the metric components are C 1 functions in these 
coordinates. 

Theorem 3.22. Assume g = g L Ug R is smooth on either side o/E ; and Lipschitz 
continuous across E in some fixed coordinate system at each point in E. Then for 
each point pinU, [JC] =0, if and only if, there exists another coordinate system, 
defined in a neighborhood of p, such that the metric components are C 1 ' 1 functions 
of these coordinates, and are related to the original coordinates by a C 1,1 coordinate 
transformation. 
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Proof. The initial coordinate system is related to an MGS coordinate system by 
a C ' transformation. By corollary |3.8| , if [JC] = 0, then the metric components in the 
MGS coordinate system are C 1 . Furthermore, since g is C 2 on each side of E it follows 
that g a f3,nv is bounded for each a, /3, 7 = 0, . . . , n — 1, therefore each components is 

Conversely, if the metric components g a p in the original coordinates are equivalent 
to a C 1 ' 1 metric in another coordinate system by a C 1 ' 1 transformation, then in an 
MGS coordinate system, which is related to the latter coordinates by a C 1 ' trans- 
formation, the metric components will be C 1 ' 1 functions of these latter coordinates. 
Therefore, by corollary |3.8|, we have [JC] = 0, and the theorem is proved. □ 



Existence of Local Lorentzian Frame Via a C > Transformation. 

Theorem 3.23. Assume g = g L Ug R is smooth on either side o/E, and Lipschitz 
continuous across E in some fixed coordinate system at each point in E. The [JC] = 
at each point of E, if and only if, for each p in E, there exists a coordinate system 
that is locally Lorentzian at p, and is related to the original coordinates by a C 1 ' 
coordinate transformation. 



Proof. If [JC] 7^ at some point p in E, then, by Theorem |3.22| , there exists no 



coordinate system defined in a neighborhood of p, and related to the original coor- 
dinates by a C 1 ' 1 coordinate transformation, such that the metric components are 
C 1 functions of these coordinates. As a consequence, there cannot exist a Lorentzian 
coordinate frame containing p and related to the original coordinates by a C 1 ' 1 coor- 
dinate transformation. 

Conversely, suppose [JC] = at every point in E. Denote the original n dimensional 
coordinate system, at a point p in E, by ). Then, by lemma |3.3|, these 



coordinates are related to a MGS coordinate system by a C 1,1 transformation. By 
Lemma |3.5| we know that gu > 0, for i = 2, ... ,n — 1, in any MGS coordinate 



system. Thus, we can choose gu = 1 and guj = at p for i, j = 2, . . . ,n — 1. 
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Now, we construct a coordinate system, denoted by u = (u°, . . . , w n_1 ), just as we 



did in Remark |3.6| , so that the metric is of the form g = diag(— 1, 1, . . . ,1). By the 
construction of this coordinates system (see the proof of Lemma |3.5| and Remark |3.6|) , 
the transformation from our MGS coordinate system to u coordinates is in the class of 
C 1 ' 1 functions. Therefore, it remains to show that g aa ,i3 — for all a, (3 — 0, . . . , n— 1. 
For i, j = 2, ... ,n — 1, we already have gij >a = where a = 1, ... ,n — 1, and 
Kij = gijfl = 0. Also, we have g 00 ^ = for (3 = 0, . . . , n — 1. 

To show gii,/? = we first consider the metric components in MGS coordinates 
denoted by g^p 5 . In MGS coordinates, gff GS = throughout E, so g^[ G a s = for 



1, . . . ,7i — l, and /Cn = gi( GS = 0. The metric component gn in w-coordinates, 



since x and u coordinates describe the same neighborhood of p, can be written in 
terms of metric components in MGS coordinates as 

. _ „MGS ^ „MGS , 1 „MGS 

9n — 9oo 9oi + ^ 9n ■ 

Therefore, since g^ G ^ = and g^jf — at p, then g u i3 = for all (3 = 0, . . . , n — 
1. □ 



Justification of the Main Theorem, Theorem |3.9| . The justification of The 



orem |3J] comes from the statement of Theorems |3.17| , |3.21| , |3.22| , and p.23| . 



3.4. The Spherically Symmetric Case in Four Dimensional Spacetime 

We conclude this chapter with another central result which considers the case of 
matching two spherically symmetric metrics across a hypersurface in four dimensional 
spacetime. The theorem below will show that the weak form of conservation across 
the shock surface, [G a/3 ~\ N a Np = 0, is implied by only a single condition when the 
areas of the spheres of symmetry match smoothly at the shock surface and change 
monotonically as the shock moves transversely to the areas of the spheres of symmetry. 
The argument in the proof is formulated so that this implication holds even when 
the shock surface is null. Smoller and Temple have already proved this theorem for 
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the non-null case in [14], and the proof, of the non-null case, can also be found in 



and ||16|| . Therefore, we will prove the theorem here only in the case when n is null. 
Our argument differs only slightly from Smoller and Temple's in that we do not make 
use of the Israel condition 

[Gqq] = [trace(fT 2 ) - (traced) 2 ] , 

but instead compute [G o] directly. However, we could prove the non-null case us- 
ing the same argument below, and replacing N with n and MGS coordinates with 
Gaussian normal coordinates. 

Theorem 3.24. Assume that g and g are two spherically symmetric metrics that 
match Lipschitz continuously across a three-dimensional shock surface S to form the 
matched metric g U g. That is, assume that g and g are Lorentzian metrics given by 

(3.66) ds 2 = -a(t, r)dt 2 + b(t, r)dr 2 + c(t, r)dQ 2 
and 

(3.67) ds 2 = -a(t, f)d? + b(i, f)df 2 + c(t, f)dQ 2 , 

where dfl 2 = d6 2 + sin 2 9dip 2 is the standard metric on the unit 2 sphere. Assume 
that there exists a smooth coordinate transformation \& : (t,r) — > (t,f), defined in a 
neighborhood of the shock surface £ given by r = r(t), such that the metrics agree on 
S. (We implicitly assume that 9 and tp are continuous across the surface.) Assume 
that 

(3.68) c{t,r) = c{i,f), 

in an open neighborhood of the shock surface S ; so that, in particular, the area of 
the two-spheres of symmetry in the barred and unbarred metrics agree on the shock 
surface. Assume also that the shock surface r = r(t) in unbarred coordinates is mapped 
to the surface f = f(t) by (i,f(t)) = ^(t,r(t)). Let N(c) denote the derivative of the 
function c in the direction of the vector N ; and assume that N(c) ^ 0. 
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Then the following are equivalent to the statement that the components of the 
metric g U g in any MGS coordinate system are C ' functions of these coordinates 
across the surface S : 

(3.69) [G aP }N a = 

(3.70) [G aP }N a Nf = 0, 
and 

(3.71) [K] = 0. 



Proof. From Theorem ^9 we have that [K\ = at each point of £, if and only 
if, for each point p in E there exists a C 1 ' 1 coordinate transformation defined in a 
neighborhood of p, such that, in the new coordinates the metric components are C 1,1 
functions of these coordinates. 



The functions c(t,r) and c(t,r), which are defined in equations (|3.66|) and (|3.67| ) 



transform as functions under (t, ^-transformations. By equation (|3.68|) we have 522 = 
c = c = g 2 2 on E, and therefore #22,0 = N(c) = N(c) = 522,0 ^ and 5-33,0 = 
N(c) sin 2 9 = N(c) sin 2 9 = g 33fi / on E. Furthermore, 522,00 = N'(c) = N'(c) = 
^22,00 an d ^.oo = N^^sin 2 ^ = N'(c)sin 2 # = 533,00 on Then we also have from 
Theorem |315] that since [JC] = and [JCtt,o] — [9tt,oo] = for t = 2, 3, it follows that 



[G aP \N a = 0. Furthermore, if [G aP ]N a = 0, then [G aP ]N a N^ = 0. Therefore, to 
complete the proof, it remains to show that [G a/3 ]N a N 13 = implies that [/C] = 0. 

To this end, we choose a smooth coordinate system (w 1 , w 2 , w 3 ) = (z 1 , 8, if) where 
d/dw 1 = n. Now extend the above coordinates to a MGS coordinate system w = 
(w° , w 1 , w 2 , w 3 ) . Notice that by the construction of a MGS coordinate system, N is 
perpendicular to the coinciding 2-spheres of symmetry, that is, N depends only on 
the time and radial components of any coordinate system. Furthermore, in MGS 
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coordinates, the metric g U g is diagonal except for (<? U #)oi = rj which is a constant. 



Since, by Lemma pT7|, fC a b = —(l/2)(g{Jg) a b,o, we have that K a b is diagonal. Therefore, 
the only non-zero components of K, are 



(3.72) 



(3.73) 

and 

(3.74) 

Thus, we have 

(3.75) 

and 

(3.76) 



£■22 = --{9^9)22,0, 



£33 = —ji.9 U 5 , ) 3 3,o- 



[/c 22 ] = 0, 



[£33] = 0, 



across E. Thus it only remains to show that [/Cn] = 0. To do this we must use our 
condition [G aP ]N a N p = 
In MGS coordinates this is 

(3.77) [G^] N a N? = [Goo] = go* [G a ] = goo [Gg] + g i [G l ] = + rj [G ] = 0. 



By equation ( |3.42|) , 
(3.78) 



[g°o] = - e Kl] = - k 2 ] - - M 



and by equation (|3.27|) we find that 



[i&J = - [lC\K tt ] - VlK\ = -g 11 [Ku + I* ) [/C 
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for t — 2, 3. Also, from equation ( |3.27| ) it follows [-R323] = 0. Now since, for t = 2,3, 

r« = -^g 08 {—gu,s + 9st,t + gts,t} = —-^g m gtt,i = —-^gu,i, 

where we have used that g 01 = I/77 and goi — V are the only non-zero off diagonal 
components of g and g~ l , and so we can say 



(3.79) 



[R l m] = -g 11 Uu - ^gu,x) [JC U ] 



Therefore, by equations (|3.78| ) and ( |3.79| ), and also using g 11 = —\/rf, g 2 2,i = n(c) 
and g 33j i = n(c) sin 2 6, we can write 

(3.80) 

[G° ] = -g 22 [Rlu] - g 33 [RU 

= V - i, 922 ' 1 ) [/Cn] + ^ - ^ 33 - 1 ) [/Cn] 



— ( - -n(r) ) [!C U ] ! J 2 ( N(c)sin 2 0-^n(c)sin 2 0) [/C 



1 

2r/~ 



cr/ 

2 r N(c)-±n(c))[K n 



csin 6*7/ 2 



1 

277" 



11 



Using equation (|3.43| 



(3.81) [Gl] = Will = + = ^ [^»] + ^ ^03] , 



t^0,l 



where we have used g^ = for i 7^ j and i, j = 2, 3. By equation (|3.28|) , for t = 2, 3, 
(3-82) [R] ot } = - [RU = 1 [i^n] + ^,1) • 



* [ G o] = ^ (N(c) + in(c) ) [/C nJ 



Therefore, 
(3.83) 

Using equations ( |3.80|) and (|3.83| ) we can write equation ( |3.77[) as 
(3.84) 



[G aP \ N a N p = [JC n ] = 0. 



cr}'- 
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Since we assumed that N(c) 7^ 0, we must have that [IC n ] = 0, and the theorem is 
proved. □ 
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CHAPTER 4 



An Exact Lightlike Shock- Wave Solution of the Einstein 

Equations 

The goal of this chapter is to construct exact, spherically symmetric, lightlike 
shock-wave solutions of the Einstein equations. We will follow Smoller and Temple's 
work from [ |15|j in which they matched a Friedman- Robert son- Walker (FRW) metric 
to a Tolman-Oppenheimer-Volkoff (TOV) metric across a timelike hypersurface. We 
will use their same matching technique here to show the existence of lightlike shock- 
wave solutions of the Einstein equations. 

We are modelling spherically symmetric lightlike shock-wave expanding into an 
static spacetime. The geometry of the region behind the shock will be that of our 
universe given by the FRW metric, and the region in front of the shock will be static 
spacetime whose geometry is given by the TOV metric. 

In this chapter we first derive the exact FRW and TOV type solutions of the 
Einstein equations, and then match these solutions along the 2-spheres of symmetry 
to get an exact, lightlike, shock-wave solution. As stated in the introduction we will 



follow [15 1, but will use the prescription given in the previous chapter to obtain our 



lightlike shock solution. 



4.1. The Metrics 

As we stated in chapter [|, physical principles determine our spacetime metric. 
The case is no different here. Both the FRW and TOV metrics are derived from 
assumptions which seem to concur with observational data. 
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The Friedman-Roberson- Walker (FRW) Metric. The FRW metric is gen- 
erally regarded as the metric which models our universe on a large scale, around 



10 s to 10 9 light years [20]. It is derived from the cosmological principle which are 
assumptions that our universe is both homogeneous and isotropic at each point. Ho- 
mogeneous means that, on a large enough scale, in any given instant of time, each 
point of space looks like any other. Isotropic means that there are no preferred di- 
rections in space; observations do not depend upon which direction we look when we 
consider these large scales. We will not derive metric here, but the derivation can 
be found in most texts on the subject of general relativity. See j5), §|, ||12|| , |18 



and |20[| just to name a few. The FRW metric is defined as 



(4.1) ds 2 = -dt 2 + -^- dr 2 + r 2 R 2 (t)dQ 2 , 

1 — kr l 

where dfl 2 = d6 2 + sin 2 6 d(p 2 denotes the standard metric on the unit 2-sphere. The 
constant k can be chosen to be +1, —1, or each giving the spatial geometry of the 3- 
sphere, flat space, and hyperbolic space respectively. The function R(t) is sometimes 
referred to as the "cosmological scale factor." 

This FRW metric was first derived by H. P. Robertson and A. G. Walker in the 
1930's, but it was incomplete in that it did not give a prediction for the function 
R(t). In 1922 Alexandre Friedmann had made some assumptions about the material 
content of the universe, and if one derives FRW metric as a solution of the Einstein 
equations with these assumptions, then R(t) can be computed (20 . 



The Tolman-Oppenheimer-Volkov (TOV) Metric. The TOV metric de- 
scribes the geometry in front of our shock surface which models a static and isotropic 



gravitational field. An especially good derivation of this metric can be found in [20 
chapter 11]. The TOV metric is given by 



(4.2) 



ds 2 = -B(f)dt 2 + A(f)- 1 df 2 + fW. 
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The TOV metric given in barred coordinates to distinguish it from the unbarred 
coordinates of the FRW metric. 



4.2. The FRW and TOV Solutions 

In this section we derive the exact FRW and TOV type solutions of the Einstein 
equations which will encompass the regions of spacetime behind and in front of the 
shock surface. 

We assume matter in spacetime is modelled by a perfect fluid that is comoving 



relative to the coordinates, that is, in free fall relative to the coordinates [20, section 



11.8]. Then the four-velocity of the fluid is given by 



(4.3) u u = v/^oo, and u l = 0, i = 1, 2, 3. 



Recall that (see equation (|2.44j )) for a perfect fluid the stress-energy tensor becomes 



(4.4) T Q/3 = pg a ? + (p + p)u a vP, a,p = 0,...,3. 

The FRW Solution. The FRW solution of the Einstein equations will be the 
solution behind the shock. Substituting the FRW metric ([4. 1|) into the Einstein field 
equations 

G = kT, 

where T is given by (f4.4|), and assuming an equation of state of the form p = p(p) 
yields the following pair of differential equations 

(4.5) P= ~ p -^ 



and 

(4.6) H 2 + k = ^^pR 



2 , ,„ 87r £ D 2 



from which we can solve for the unknown functions R(t), and p(t) |15| , [16 
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As a simplification we restrict ourselves to the case k = so that the FRW metric 
in (|4.1p becomes 

ds 2 = -dt 2 + R 2 (t)dr 2 + r 2 R 2 (t)dn 2 , 
and the differential equations (|4.5| ) and Q4.6Q become 

(4.7) 

and 

(4.8) R 2 = ^-pR 2 . 



(4-9) R = ±( 8 -^) R, 



Rewriting the equation ( |4.8|) as 

1 2 

3 ) 

and then substituting into equation ( |4.7| ) we get 

P 



(4.10) P = ~PT 



Notice that the ± signs in equation ( |4.9j ) correspond directly to the =F signs in equa- 
tion ( [4.101) . We can solve ( |4.10|) explicitly when the equation of state p = p(p) is 
given. Indeed, solving ( [4.10| ) for dt yields 

(4.11) dt = T~, d P , 

which we integrate to get 

di 



(4.12) t-t = T 

Then using by equation ( f4.9|) , and 



l>0 



• dp dR , , r— — —dR 
(4.13) fl = £_ = T( , + p) ^i^_ 

we have 

dR —dp 



(4.14) 



R 3(p + pY 
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(i<) 



It follows that equation ( |4.14| ) has solution 



(4.15) 



R = Rq exp 



v, 3(€+p(0) 

The TOV Solution. Now we give the the TOV solution which encompasses the 
region in front of the shock, and represents a general relativistic version of static, 



singular isothermal sphere |15 



We proceed in a similar manner at the FRW metric by substituting the TOV 
metric ( f4.2|) into the Einstein equations 

G = kT, 



where T is given by ( |4.4| ). In this case, the substitution (see pOfl ) gives the following 



(4.16) A(f) = ( 1 - 



(4.17) — -=4vrf 2 p, 

dr 

and 

Here the unknown functions, p(f), p(f), and M(f), depend only on f. Analogous to 
the FRW solution we assume an equation of state p = p(p) for the TOV metric so 
that the differential equations ( [4.17| ) and ( |4 . 1 8| ) give solutions for M(f) and p(f). The 



function M(r) denotes the mass inside radius r, and can be written as 
(4.19) M(f)= / 4vr£ 2 p(0^- 

To find the metric component B{f) we look to the equation for hydrostatic equilib- 
rium ||20| , equation 11.1.8], 



B'ir) p'(r) 
4.20 = -2^-. 

v ' B p+p 
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Now we further restrict the equation of state for the TOV metric to be of the form 
(4.21) p = ap, 

where a is a constant, and also assume the energy density is given by 
( 4 - 22 ) P = ^ 



with 7 a constant. With these assumptions equation (|4.19 ) becomes 



(4.23) M(f) = / 4< 2 p(0 d£ = 47T7f 



r 

2- 



Combining equations Q4.21J) - (|4.23|) with ( |4. 18|) gives 



(y4 ' 24 " > 7 2vr£ (l + 6a + a 2 

and inserting equation (|4.23| ) into equation ( t4.16|) yields 



(4.25) A = 1 - 87r£ 7 . 



Now, we find B by inserting equations ( [4.21 ) and ( f4.22| ) into equation ( [4.20 ), which 
becomes 



(4.26) 

Then solving we get 
(4.27) 



dB 2a dp 

~B = ~(1 + 



\ 4ct/(1+ct) 



4.3. Matching the FRW metric to the TOV metric 

In ||14|| Smoller and Temple derive a coordinate transformation that takes (t, f) — > 
(t, r) so that the FRW metric ( |4.1| ) and the TOV metric (|4.2|) match Lipschitz con- 
tinuously across a shock surface S. The same derivation can also be found in Q 
and ||16|j . Here we summarize their procedure here as briefly as possible highlighting 
their results which we will use to construct a lightlike shock-wave solution. 
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Since our goal is to be able to apply Theorem |3.24| we begin by letting 



(4.28) 



which ensures that 



r(t, r) = R(t)r, 



f 2 dQ 2 = R 2 r 2 dQ 2 . 
A consequence of equation ( |4.28| ) is that we can write the following 



(4.29) 
(4.30) 
(4.31) 



dr = Rdr + Rrdt, 



1 ,- R , 
dr = —dr — —rdt, 

R R 

f Rr 



R R 

Now, it can be shown that in (t, f)-coordinates the FRW metric ( |4.1|) is written as 



(4.32) ds 2 



.#2(1 _ ^ pR 2 r 2 ) dt 2 + R 2 df 2 - IRRfdt dr \ + f 2 dQ 2 



R 2 - kf 2 { " v " 3 
At this point a mapping t = t(i, f) is constructed to eliminate the cross term dt dr in 



equation ( |4.32| ). This is done first for a general metric given by 
(4.33) ds 2 = -C(t, f)dt 2 + D(t, f)df 2 + 2E(t, f)dtdf, 



and choosing if) = ip{tj r ) such that 
(4.34) 
so that 



(4.35) 



dt = if)(t, r){C(t, f)dt - E(t, f)dr}, 



is an exact differential. This means that equation (|4.33| ) can be written as 

(4.36) ds 2 = -(^C-^dt 2 +( D + ^P) df2 - 
Therefore, the FRW metric in (i, f)-coordinates is 

(4.37) ds 2 = R2 l _ k _ 2 {-(^cy'dt 2 +(d + ^\ dr 2 \ + f 2 dn 2 , 
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where 

(4.38) C = R 2 {1 - ^-pR 2 r 2 }, D = R 2 , and E = —RRf. 

Now we can finally define the shock surface at which the FRW metric ( |4.1| ) and 
the TOV metric ( |4.2|) match Lipschitz continuously. Indeed, using equation ( (4.38[) 
write 

, \ r. E 2 r.2 R 2 R 2 f 2 n2 R 2 R 2 r 2 

( 4 - 39 ) D + — = R 2 + - , 8 , . =R 2 + 



C i? 2 (1 - ^f-pf 2 ) l-*fpR2 r 2- 

and equate the df 2 components in the TOV metric ( |4.2|) and the FRW metric in 
(i, f)-coordinates (|4.37). Then it follows that 



(4.40) 



- ^ P R 2 r 21 
and using the differential equation ( |4.6| ) we get 

(4.41) M(f) = yp(t)f 3 , 

which implicitly defines the shock surface E. To express S in (t, r)-coordinates we 
use the transformation f = R(t)r. For the shock surface S to remain in the domain of 
definition of the FRW metric we must assume that 1 — kr 2 > for k > 111 411. We also 



must obtain the conditions under which the function ip defined in equation ( 4~3" 
which determines t from (t, f), can be found uniquely. Since the dt 2 terms must match 
on the shock surface, it follows that 



(4 - 42) wh^ik = B(f) 

must hold on the shock surface. Then tp can be determined by the partial differential 
equation 



(4.43) 



cvv- + i#t = /(t,f,v>), 
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14|, and [16 



where C(t, f) and E(t, f) are given in equation ( |4.38| ). As shown in [[ 
the partial differential equation ( |4.43|) with initial value given by ( |4.42| ) can be solved 
uniquely in neighborhood of a point on the shock surface provided 



(4.44) 



df C 
~dt ~E' 



Here r denotes the speed of the shock surface, and equation ( |4.44| ) is the condition 



that the shock surface be non-characteristic at a point [|14| . 

Now that we have the shock surface, we restate another result of Smoller and 
Temple, proposition 1 in | 



[ |14|| , and | |16|| , which will be used for computations 
involving the transformation from (t, r)-coordinates to (i, f)-coordinates and back. 



Proposition 4.1. On the shock surface given by equation ( |4.41| ) ; the following 
identities hold: 



(4.45) 

(4.46) 

(4.47) 

(4.48) 
(4.49) 



1 


( AE 2 \ 

B { 1+ C 2 ) = 


c = 


R 2 A, 


E 


-Rr 


c 


A ' 


/■;- 


-A+(l- kr 2 ) 


c 2 ~ 


A 2 


R 2 r 2 = 


-A + (l- kr 2 ). 



B 
A< 



1 - kr z 



In addition to using proposition |4.1| to transform quantities in (t, r)-coordinates 
to and from (i, f)-coordinates we will also need the following proposition. 



Proposition 4.2. On the shock surface given by ( |4.41| ), the following partial 
derivatives can be written as 
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Proof. From equation (|4.35|) we can write 



(4.52) 



from which, follows 



dt = (ipC)- 1 dt+ —df, 



(4.53) 



and 



dt 
dl 



By equation ( (4.30[) we have 



dt _E 
df~C' 



1 ,_ Rrdt ,_ 



and so, using equation ( |4.47| ) 



dr _ Rrdt _ AE x 
ffi ~ lldt~ RC^ ' ' 



Lastly, using equations ( |4.28| ), ( $AT\) , ( |4.48| ), and ( |4.53| ) yields 



<9r 9 / r 
df df \R 



1 i?r dt 1 _Rr 



1 

R 



A- 



E\ E 1 



C CR 



1 

7? 



l - 



A + (1 - Ax 2 



1 — A;r 2 
RA ' 



□ 



We summarize the conditions under which we can match the FRW and TOV 
metrics, given in ( |4. 1[) and ( f4.2|) respectively, in the following theorem, which is a 
restatement of theorem 6 in [0], and theorem 5 in WM and |T6| . 



Theorem 4.3. Let the point (t ,r ) satisfy 

(4.54) M(f) = — p(t)f 3 , 

3 
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and let ( |4.54| ) define the shock surface f = f{t) in a neighborhood of (t Q ,f ). Further- 



more, assume 

(4.55) R(t)r = f, 

requiring the barred and unbarred coordinates to be equivalent on the spheres of sym- 
metry of the FRW and TO V metrics, given in ( |4.1| ) and ( |4.2| ), and the shock surface 
in (t , r) - coordinates be given by r = r(t) = f(t)/R(t). Then, when 



1 - kr 2 > 0, 
A(f ) 0, 

and 

df C Rr 
T = ~dl^E = ~~A' 

hold at t = to, we can define the coordinate i so that there exists a smooth regular 
transformation on a neighborhood o/(to,^o) that takes 

(*,*•)— >(f,f), 



and the FRW and TO V metrics, given in (4.1) and (4.2), match Lipschitz continu- 
ously across the shock surface r = r(t). 

4.4. The Conservation Condition 

The goal of this section is write down a formula that expresses the weak form 
of conservation of energy across a spherically symmetric lightlike shock surface con- 
structed in the manner of the previous section. We should note that, so far, in this 
chapter we have merely recapitulated results from H, [|14||, and ||16||, that is, we have 



not stated anything that has to do with our shock surface being lightlike or not. It 
is in this part of this chapter that our results begin to diverge from that of Smoller 
and Temple in the sense that we are dealing with the lightlike case. However, we are 
still using their work as a model to obtain our own results. 
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The Lightlike Surface. Assume that equation ( f4.41| ), 

47T 



(4.56) M(f) = ^-p(t)f 3 , 
defines the surface given by the level curve 

(4.57) <p(t,r)=r-r(t) = 

in a neighborhood of (to^o) which satisfies equation (4.56). We can also write (p in 



(t, r)-coordinates as 

(4-58) <p(t, f) = J - r (t(t, f )) . 

Then we can compute the normal to (p — by 

(4.59) dip = n dt + n\ dr = —rdt + dr 
to find 

(4.60) n = — r, and n\ = 1. 

Although, the computations are a little more complicated, from p> in equation ( |4. 58 ) 
we can write 



- T — tRt 

dV = n»dt + n l df=-— ) dt + —df 



which yields 



T 1 tRt 

(461) ^ = -Wcy and " 1 = ^ + ^i- 

We are assuming our surface is light like, hence 

1 — kr 2 



(4.62) (n, n) = g^np = g 00 (n Q ) 2 + g 11 ^) 2 = -r 2 + 



R 2 

which implies that 

1 — kr 2 



(4.63) r 



2 



R 2 

Notice that if we compute 

(n,n)=g 00 (n ) 2 + g u (n 1 ) 2 
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using equation ( |4.61| ), and the transform the result into (t, r)-coordinates we get 



equation (|4.62| ) 



The Transverse Vector. Equation ( [4.62| ) implies that n lies in that tangent 
space of the surface, and therefore we must choose a transverse vector N which 
satisfies equations 



(|3.5| ), and (|3.7f ). Now for any N we have 
(4.64) (N, n) = g 00 N n + g n N l n l = N°n + N l rn = -fN° + N l = r] ^ 0, 



In light of equation ( |4.64| ) , we choose N so that 

(4.65) iV° = 0, and N 1 = rj ^ 0, 



or, lowering indices, 



(4.66) 



N = 0, and iVi = gnN 1 = r\ 



R 2 



1 — kr 2 



Lemma 4.4. The transverse vector N defined in equation ( [4.65 ) satisfies 



(4.67) 



[<N,X a >]=0, 



where {X a }^ =1 is a basis for the tangent of the surface at (to, ro), and 



(4.68) 



(N,n) = (N,n)=r^0, 



where (N, n) is written in (t,r)- coordinates. Furthermore, we have that 



(4.69) 



[<N,N>] = 0. 



Proof. To show equation (|4.68|) we first transform Nq to iVo, and Ni to Ni. 



Using propositions [4.1| and |4.2| we have 



(4 . 70 ) ^ = ^1 + ^1 = ^^W- = ^(W, 



and 
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Then, using equations (|4~2]), ( |L45[ ), ( |L47[ ), and ( |L61| ) we find 



(N,h)=f°N n + g n N 1 n 1 



fRr 



(4.72) 



B- l Ar E B(l-kr 2 ) 



»? I — + 1 + 



.4 



.4 



.4 



Then equations ((164]), (f!6|),and (gTFg ) prove ( gjgD that holds. 

For equation (|4.67f) we must choose a basis for the tangent of the surface at (t , fo) 

Let 



Xi = n, X 2 = ^, andX 3 = ^ 



Then by equations ( |4.64|) , and ( |4.72|) we have 



[(N,X 1 )] = (N,n)-<N,n>= V ~ V = 0. 



Now, 



[<N,JQ]=0 



for a = 2, 3 since 
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To show equation Q4-69Q holds we compute (N, N) in (t, r)-coordinates. From 
equations (|4.70|) and ( |4.71| ), and the relations in proposition [4.1| we see that 

(N,N) =g 00 N N + g 11 N 1 N 1 



(4.73) 



1 R 2 A 2 



5(1 -Ax 2 ) 2 
R 2 



-A+l-kr 2 )B(l-kr 2 ) R 2 

T 2 A + T 



R*_ R*_ 

1 - kr 2 A A 

R 2 



1 — kr 2 



We also have 

1 _ kr 2 / R 2 \ 2 R 2 
(4.74) {N,N) = g 00 N N + g 11 N 1 N 1 = ^^ V -^—) = V 2 



R 2 



1 — kr z 



1 — kr 2 



and thus 



[(N,N)] = (N,N>-<N,N>=77 2 



R 2 



1 — kr 2 



- rj 



R 2 



1 — kr 2 



0. 



□ 



The Conservation Condition. The conservation of energy across the surface 
defined by equation ([4.41| ) is given by 



(4.75) 



[T^N a N p ] = 0, 



where T a/3 is the stress energy tensor for a perfect fluid given in equation ( |4.4| ). 
For comoving coordinates, see equation (fO|), the jump condition ( |4.75| ), using equa- 
tions (|4.66|) and ( |4.70|) , becomes 

B 



(4.76) 



R 2 , ^ O _T / 1 1 



[T^N a N p ] = (p - p)|N| 2 + (p + p)N 2 -(p + p)^ 

= {P-P)v 2 -, - '/' - /';'/ < 

1 — kr A 

it ■■■ n n + o \ 



A \-kr' 



2 p2 / P + P P + P 



80 4. AN EXACT LIGHTLIKE SHOCK- WAVE SOLUTION OF THE EINSTEIN EQUATIONS 

where we have used 



N 2 = ( VR 



A E 



(W 1 



1 - kr 2 C 
R 2 A 2 (-A + l-kr 2 )B(l-kr 2 ) 



[1 - kr 2 ) 2 



7] 2 BR 2 



A 2 
1 



A 



1-kr 2 A 



When the jump condition ( |4.76|) holds, it follows that on solutions of the Einstein 
equations G = kT, this is equivalent to the jump condition on the Einstein tensor 



(4.77) 



[G af} N a N p ] = 0. 



Therefore, all the equivalencies in Theorem |3.24| hold which implies all of the equiva- 
lencies hold in Theorem [3.9| . Thus, there exists a lightlike shock-wave solution of the 
Einstein equations when Theorem [O] and equation ( |4.76| ) holds. In the next section 
we give an exact such solution. 



4.5. An Exact Lightlike Shock- Wave Solution of the Einstein Equations 



The jump Conditions. We begin by assuming the suppositions of Theorem 4.3 



are satisfied, the jump condition ( |4.76[ ) holds, k — 0, and the equation of state for 
the TOV metric is given by 

P = 

for some constant a. Furthermore, also suppose that the TOV solution is given 
by (|4.22 ) - Q4.25 ), and ( 4.27 ). On the shock surface given by 

An 

M{f) = — P {t)f z 

we solve for p via the coordinate transformation f(t) = r(t)R(t) to find 



(4.78) 



3 M 



37 



47rf(t) 3 f(ty 



3p. 



4.5. AN EXACT LIGHTLIKE SHOCK- WAVE SOLUTION OF THE EINSTEIN EQUATIONS 81 




Figure 1. A graph of a = H(a) defined in equation (|4.80|) . Note that 
H(l) = 1 and H(0) = 0. 

Now, we are able to compute the pressure p of the FRW metric. Indeed, substituting 
k = 0, p = op and p = 3p into the jump condition ( |4.76|) yields 



(4.79) 



where 



P = °P, 



1 (a + l-A 
a= 3 [— 4— 



3(a + l)' 



Analogous to [14 we can state this relation as 



(4.80) 



a 



a'- 



18a + 25 + -a- - = H(a) 



2 2 2 

As in non-lightlike case, see | |1 5| | , within the region < a, a < 1, we have H(0) = 0, 
H'(a) > 0, and a < a; see figure [1]. However, in our case H(l) = 1, which is not 
true in ||15|| . As we shall see below, having to satisfy the lightlike condition ( [4.63D 
will restrict cr to a single value. So far, we have that when H(a) = a, the jump 
condition (|4.76|) holds, and so we can apply Theorem |3.24| which means that we also 
get the equivalencies of Theorem |3.9| . This means we can say that the matched FRW 
and TOV metrics form a lightlike solution of the Einstein equations if the shock speed 



in (t, r)-coordinates satisfy the lightlike condition (|4.63). Therefore, we must compute 
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the formulas for p(t), R(t), the shock positions r(t) and r(t), the shock speeds f(t) 
and r(t), and then find the appropriate values of a so that equation ( |4.63| ) holds. 

The Shock Solution. We begin by substituting p = ap into the FRW solu- 
tions (pLTll) , and (gig) to get 



(4-81) dt = T^=±-—p- 3/2 dp, 

and 

dR 1 dp 

4.82 — = ; - — . 

y J R 3(1 + a) p 

From fl47p we can write r in terms of p, 



(4.83) f = ^p- 1/2 , 
and then differentiating to get 

(4.84) p~ 3/2 dp = ~^^ df - 



Notice that we can directly substitute this expression for p 3 ^ 2 dp into equation ( |4.81| ) 
which yields 

(4.85) dt = =F 7- 1 r ^ — ~ dv . 
Now integrating gives the shock position 

(4.86) f{t) = ±Vl87rS 7 (l + *){t- *o) + r , 



and in conjunction with (|4.78| ) we can compute the FRW energy density as a function 
of t, 

(4.87) p( f ) 3 ^ 3 ~> 



f(t) 2 (±^T8i37(l+ff)(f-f„)+f„] 2 ' 
Now, we can solve the differential equation (|4.82j) for R(t) to find 

✓ v -1/3(1+*) , -/.xv 2/3(1+<t) 

(4.88) =*(-^) • 
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from which we obtain 

/f(f\ \ -2/3(l+<0 

r(t) =f{t)R(t)- 1 = f{t)Ro-l 1 

(4.89) 



r 

(1+3ct)/(3+3ct) 

f R -l 



r J 

The Shock Speeds. To compute the shock speeds in (t, f )- and (t, r)-coordinates 



we differentiate equations (|4.86|) and ( [4.89|) with respect to t, and get 



and 



1 + 3a / a 

(4.91) r 



R(t) V l + 6a + a 2 ' 
The Lightlike Condition. In Smoller and Temple show that the shock 
speed, relative to the FRW particles, is given by 



(4.92) s{a) = (1 + 3a) 



1 + 6a + a 2 ' 

where we use a = H(a), which is defined in equation ( |4.80| ), to make this expression 



a function of only a. The function s(a) is valid in our case as well, where the only 
deviation from Smoller and Temple's argument, given in [I5L lemma 1], is that we 



must use our version of H{a) which is slightly different than the one they derived. 

For our lightlike solution to be valid we must satisfy the extra condition which is 
the lightlike condition given in ( f4.63| ) where k = 0. Therefore, we must have that 



(4.93) s(a) = 1. 
Solving s(a) = 1 for a, see figure |2], gives 

(4.94) a = a 2 « 0.63442. 



We can now characterize our lightlike shock-wave solution which we do in the 
following theorem. 
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l-s(a) 




±0.4- 



±0.6- 



FlGURE 2. A graph of l — s(o~) defined in equation ( |4.94j ). The solution 
of 1 — s(er) = gives a value of a = cr 2 ~ 0.63442. 

Theorem 4.5. Assume an equation of state of the form p = p(p) for the TOV 
metric, and p = op for the FRW metric. Also, assume that relation between a and 
a, given by 

1 



a = -V9a 2 - 18(7 + 25 + -a - - = H(a) 
2 2 2 v 1 



holds, where we have taken k = 0. Then the TOV solution given by 



p = 7/r 2 



M{r) = 47T7r, 
A = 1 - 8tt£7, 



and 



B = B n {- 

r 



- \ 4a 2 /(l+a 2 ) 



will match the FRW solution given by 



P(t) 



37 



and 



[±v / T8^(l + a 2 )(t - t ) + f ] 2 ' 

«( () ^(f) 2/3(1+ " 
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across the shock surface 



f(t) = ±v/l87T07(l + a 2 ){t - t ) + f , 

such that conservation of energy and momentum hold across the surface. We have 
used o~2, and H{a 2 ) = o 2 , where o 2 is the solution of s(o~) = 1, and s(a) denotes the 



shock speed given in (|4.92|) 



The Lax Shock Conditions. The final question we would like to answer is how 
do we classify the shock solution given in Theorem [4.5| with respect to the shock 



conditions first given by Smoller and Temple in ||15|| ? As in ||15|| we only consider the 
case when the pressure and density are greater behind the shock wave. Since p = 3p 
this means that the FRW region is behind the shock, the TOV region is in front of the 
shock. Therefore, for our solution we take the plus sign for expression of R in 
and the corresponding signs in equations (|4.10| ), ( |4.11]) , and (|4. 12 ) . 



The shock speed given in (|4.92|) is coordinate dependent, and is computed in 



a locally Minkowskian frame |15| , |16|| . A coordinate system is called locally 
Minkowskian at a point p if g a p{p) = diag(— 1, 1, 1, 1), but not necessarily Lorentzian 
where it is also required that g a p n {p) = 0. Also, note that since we are working with 
the radial component exclusively we only need to consider a locally Minkowskian 
frame in the (t, r)-coordinates. We denote our locally Minkowskian coordinates by 
(t, f), and they can be obtained by letting r = ip{f) and choosing ip so that ip' = l/R 2 , 
which implies 

ds 2 = -dt 2 + R 2 dr 2 — > ds 2 = -dt 2 + df 2 . 

The shock speed s(a) is determined in a locally Minkowskian frame comoving with 
the FRW metric. 

Recall that we have chosen a = a 2 ~ 0.63442 to satisfy f = l/R, and since 

dr 1 df 
dt R dt' 
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it follows that 

s(o- 2 ) = 1, 

in (t, f)-coordinates. The characteristic speeds behind the shock are equal to the 
sound speed ±y/a in the (t, f)-coordinate system, since the FRW fluid is comoving 
with respect to (t, r )-coordinates ||15|| . That is, the speeds of the characteristics 
relative to the FRW fluid are given by 

*frw = = 

In our case, with the shock moving outward with respect to r and f, the Lax char- 
acteristic condition, which says that the characteristic curves of the shock family 
impinge on the shock from both sides and all other characteristics cross the shock, 
hold if and only if 



(4.95) A+ ov < s < X 



FRW 



where X^ov denotes the corresponding characteristic speed on the TOV side, or the 
front side, of the shock. Then in our lightlike case we have, 



*frw = V^~ V0.63442 = .79650 < 1 = s(a 2 ), 

which means that our shock solution does not satisfy the Lax characteristic condition. 
We can also conclude that 

(4-96) ~X FRW < X+ RW < s(a 2 ). 



Now we would like to know how 3(02) relates to the characteristics X TOV . In [15 
Smoller and Temple show that 



(4-97) X+ ov (a) = - je== —V^. 

ya z + 6<T + 1 — la 

Then for a = a 2 ~ 0.63442 we have 

Xrovfa) « -0-45040 < 1 = s(a 2 ), 
see figure |^. Therefore, we have that 
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1- 

V 

A, T0V 
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and 



Figure 3. This graph of Aj y(cr), defined in equation (|4.97|) , shows 
that Aj oy (o- 2 ) < < s(er 2 ) = 1, where a 2 « 0.63442. 



A Toy < \ TOV < s(a 2 ). 



Thus both sets of characteristics cross the shock since its speed is greater than the 
characteristic speeds on each side of the shock. 

4.6. Comparison with the Smoller-Temple Sublight Shocks 



In [ |15| Smoller and Temple show that in the limit as their subluminal shock 
solutions tend to the speed of light, 



a 



a 2 = 0.745 



However, we have shown that an actual shock solution moving at the speed of light 
the value a 2 = 0.63442. This yields the unexpected result that the solution is not 
equal to the limit of Smoller-Temple subluminous solutions as they tend to the speed 
of light. 
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CHAPTER 5 



Conclusions and Summary of Contributions 



We have generalized the work of Smoller and Temple given in [|14|] and |15| in the 
sense that we can include lightlike shock surfaces. In defining a more general second 
fundamental form, by replacing the normal vector with a transverse vector satisfying 
the jump conditions ( |3.4|) and (|3.5| ), we were able to overcome the breakdown of 
the standard second fundamental form for lightlike hypersurfaces. Then we used 
this generalized second fundamental form, in an analogous way to how Smoller and 



Temple applied the standard second fundamental form in [|14j| , to obtain a theory for 
the lightlike case. In the process, we introduced a modified Gaussian Skew coordinate 
system. This theory yielded the unexpected result of having to include the condition 
that the metric be C 2 on the spacelike subspace of T p E. 

Then we were able to construct the an exact, spherically symmetric shock- wave 
solution of the Einstein equations which propagates at the speed of light. No quan- 
tities are moving at the speed of light except the shock. Although, our solution was 
consistent with Smoller and Temple's results in [T5|, there was the unexpected dif- 



ference in the value of a^- In the non-lightlike computation given [|15| the value for 
a in the lightlike limit was 

a 2 « 0.745, 
while our lightlike computations yielded 

a* « 0.63442. 

From this we conclude that the limit of the Smoller- Temple subluminal solution as it 
tends to the speed of light is not equal to our solution propagating at the speed of 
light. 
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Furthermore, we showed that in this exact solution the pressure and density are 
finite on each side of the shock throughout the solution, the sound speeds, on each 
side of the shock, are constant and subluminous. Moreover, the pressure and density 
are smaller at the leading edge of the shock which is consistent with the entropy 
conditions in classical gas dynamics [|7], |13|] . 



5.1. Summary of Contributions 

In the lightlike shock matching theory of Chapter |3] we applied the notion of a 
generalized second fundamental form given by Barrabes and Israel in |l| into the 
shock matching framework of Smoller and Temple ||14| , |15|| . The generalized second 



fundamental form K, in n was given in a scalar form, and we modified it so that K, is 
mapping that takes tangent vectors on the surface to tangent vectors on surface. Then 
we based the analysis on a Gaussian skew type coordinate system which we had to 
modify so that we could incorporate the generalized second fundamental form into the 
shock matching theory. Most of the supporting lemmas, and each of the theorems 
in ||14|| had to be modified and proved again in the context of this new Modified 
Gaussian Skew (MGS) coordinate system and the generalized second fundamental 
form. This yielded the unexpected result of adding the extra condition in our main 
result, Theorem |3.9| , that the metric inner product on the spacelike subspace of T P E 
to be C 2 in order for the Rankine-Hugoniot jump conditions to hold. This condition 



was already satisfied for the spherically case in Theorem |3.24 



Then in Chapter [| we constructed a new exact shock-wave solution moving at 
the speed of light. In order to do this we had to construct appropriate transverse 
vectors N, and apply our extension of the subluminal theory to the specific example 
of the matched FRW/TOV metrics. This included significant modification of the 



equations expressing the Rankine-Hugoniot jump conditions given in equation ( 4.76 ). 
This yielded the difference in the limit of the Smoller- Temple subluminal solutions as 
they tended to the speed of light with the actual solution propagating at the speed 
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of light. Lastly, we incorporated the modifications of the sub luminal shock solutions 
to show that our new exact solution was a crossing shock. 
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